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APPENDIX A 


GARBLE 


A.l INTRODUCTION 
| In its simplest form radio frequency garble occurs 
whenever two signals arrive at the same receiver simultaneously 
or close enough in time that the signals interfere. When 
the condition persists with many repetitions of the overlapped 
signals, it is synchronous garble. Figure A-1 illustrates 
three overlapping replies from aircraft #2, #3, and #4 during 
one 20.3 us transponder reply period. This analysis is 


concerned only with replies on 1090-MHz since the 1030-MHz 


interrogation frequency is relatively uncluttered. The 
1090-MHz replies are omnidirectional from all aircraft in 
any radar's mainbeam (and nonsuppressed sidelobes) causing 
the BCAS to receive replies from each aircraft replying to 
each radar (passive garble). In a low density environment 
with only one or two interrogation overlaps at a maximun, 
relatively simple tracking algorithms can separate or 


"de-interleave,"' the fourteen discrete pulses of each signal 
during the 20.3 us interrogation period. But in dense 


traffic with more than 5 or 6 overlaps, passive tracking 
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becomes next to impossible. The same is true for active 


tracking at a somewhat higher level of garble. 


Active Garble 

The ground radar system becomes garbled when two air- 
craft are at the same slant range from the interrogating 
antenna and at the same azimuth angle. This is due to the 
system resolution of targets in range and the directional 


antenna resolution in azimuth. Active airborne interrogators 


normally employ nondirectional antennas due to the high cost 


and complexity of airborne directional rotating antennas or 
phased-array antennas. The active BCAS receives garbled 
replies from any target whose range from the BCAS is within 
half a message length (10.15 us or 1.67 nm) of the target of 
interest. This is a spherical shell one message length thick 


as depicted in Figure A-2. 


As DABS is deployed, the active synchronous garble 
problem is eliminated, since DABS replies are only obtained 
from those aircraft which have been discretely interrogated. 
Figure A-3 shows active BCAS synchronous garble in the LA 
Basin for the three densities at each 10 degrees around the 
target #399 for a 10-nautical-mile-range system. The total 
number of aircraft within an annulus of 10 + 1.67 nm of the 
BCAS aircraft is counted. The number of possible target 


tracks required to continue proper tracking through the 
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maximum number of overlaps was plotted for 5, 10, 15 and 20 
nm active systems. An average of these maximum overlaps was 
determined for different numbers of tracks, and this was 


plotted in Figure A-4. 


An analysis was made of the effectiveness of the 
whisper-shout technique on reducing the active garble. This 
is a multilevel (8 steps in the simulation) suppression- 
interrogation pair of pulse trains where the suppression is 
at the level of the previous interrogation in order to acquire 
only those aircraft in the next higher "shout" annular region. 
The simulation assigned each aircraft in the Basin a unique 
transponder sensitivity to approximate the actual sensitivity 
distribution determined by the MIT Lincoln Laboratory (13). 
The result of the whisper-shout simulation was that a four-to- 
one reduction of active garble level was achieved by grouping 
the aircraft into whisper shout levels according to range and 


sensitivity. Further reductions have been proposed using a 


directional antenna. 


A.3 PASSIVE GARBLE 

The passive garble region also occurs for a region in 
space where the signals from two aircraft are recevied over- 
lapped. But instead of similar slant range from the BCAS as 


with active garble, the passive garble geometry requires the 
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two aircraft to be in the same ground radar main beam and 
have equal reply times of arrival at the BCAS aircraft. 
Figure A-5 illustrates this phenomenon. The radar is at the 
far right and sends out an interrogation at time T = 0 (in 
units of message length of the reply, or 20.3 us). At T = 10 
target #1 replies, and its signal expands spherically in 
space and is illustrated at T = 12 and T= 18. Target #2 
replies to the same interrogation at T = 16, and its reply 
expands spherically and is also illustrated at T= 18. The 
intersection of these replies has been darkened in the figure 
to show the overlap region at a snapshot in time. Obviously, 
along the line from the radar through the targets beyond 
target #2, the signals remain completely overlapped. Near 
the hyperbolic boundary of the overlap region very little of 
the two signals are overlapped. It may be seen from the 
illustration in Figure A-6 that the overlap region grows in 
space as the two replies expand, and the envelope of the 
region is a hyperboloid of revolution about the radar-target 
line extension. In Figure A-5 it is shown that BCAS #1 will 
receive both replies in the clear, while BCAS #2 will receive 
a partially overlapped reply. As the spacing between targets 
#1 and #2 is varied, a family of hyperboloids is formed. 
Figure A-7 illustrates typical limits of the reply overlap 
regions, designated "garble contours." In this illustration 
target #2 is at the center of the figure, and target #1 is 


moved to the right forming hyperbolas which curve 
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progressively further to the left. Note that the vertical 
line represents a garble region (the entire left-hand half 
space) for the case where the two targets are one full message 
length apart. Also the asymptote is plotted for the case 
where target #1 is at or near the radar and #2 is at a great 
distance away. When the aircraft get closer than one message 
length (3.34 nm) the hyperboloids bend toward the radar until 
the targets are within one-half message length, at which point 
the entire space is overlapped, and the radar itself receives 
garbled replies. In this case the garble contour is a 


straight line toward the radar. 


A practical application of the difficulties involved 
in passive tracking is illustrated in Figure A-8. Three 
radars are available for tracking purposes, each at slightly 
more than 10 message lengths (say 35 nm) from the target of 
interest at the center, and at 120 degrees from each other 
with respect to the target. Assume the target is near a 
moderately busy airport with other aircraft making approaches 
and departures. If there is an additional aircraft along the 
line of sight to each radar from the target within about one 
message length from the target, then nowhere in space is the 
reply from the target seen in the clear with respect to all 
three radars. In roughly two-thirds of the airspace, one of 
the radars produces a garble of the target reply. In the 


regions adjacent to the main beam of each radar, the other two 
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radars cause the replies of the target to be garbled by air- 
craft in the main beam. And in the small area adjacent to the 
target, all three radars cause the other aircraft's replies to 


garble the target aircraft's reply. 


The foregoing discussion assumes that two aircraft 
whose replies are garbled are aligned with the radar eliciting 
the replies. This assumption is supported by the fact that 
two aircraft separated by an arbitrary distance are more 
likely to lie within the main beam at the same time, and thus 
to produce garbled replies, if they are approximately in 
alignment with the radar. Deviation of alignment of the 
aircraft with the radar causes the garble contours to become 
skewed with respect to the main beam. If the two aircraft are 
located within the main beam, at the same slant range from the 
radar, and within one message length from each other, their 
replies will be garbled everywhere. This occurs because the 
two are interrogated at the same time, thus allowing the first 
pulse of the reply sequence from one to reach the last pulse 
of the reply sequence from the other before it is transmitted. 
Although the skewing phenomenon causes some spreading of the 
distribution of garble, the assumption of target alignment 
with the radar is shown to be sound in a subsequent section 


that deals with a dynamic system simulation. 
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A.4 STATIC SIMULATION OF LA BASIN PASSIVE GARBLE 


The passive garble contours for target #399 in the 
LA Basin were plotted with respect to each of 9 radars 
(LAX, Burbank, El Toro, Long Beach, March, Norton, Ontario, 
San Pedro, Santa Ana) and are contained in Attachment I. 
Figure A-9 is the plot of hyperbolic overlap regions for the 
LAX radar (ASR-7) for the full 1982 density. Each contour line 
indicates a target in the main beam with #399 causing a region 
in space to be garbled. In the sector toward the LAX radar 
(toward the center of the axes) there are 5 overlaps due to 
targets within one half message length of #399. In the opposite 
sector there are 19 overlaps of #399. An azimuthal distribution 
of passive overlaps was derived from Figure A-9 and plotted 
alongside the corresponding active overlaps in Figure A-10. For 
comparative purposes the active BCAS garble is seen to be four 


times that of the passive. 


The LAX radar is only 20 nm from #399, and therefore 
relatively few aircraft share the main beam with #399. The 
radar interrogator at Ontario, on the other hand, causes some- 
what of a worst-case condition for garble, being some 50 nm 
away and with many more targets in the beam width capable of 
garbling aircraft #399. The symbol "R" on the right hand side 
of Figure A-11, just above the x-axis, indicates the radar 
location. The points in line with the radar and #399 are 


targets illuminated by the main beam as the radar sweeps 
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past #399 with the 4-degree beam centered on #399. The 
opposite situation is illustrated in Figure A-12; where the 
Long Beach ASR is very close (about 8 nm), the beam has not 
spread very far, and the data show no garbles on the north 
side and only three on the south side. This may be contrasted 
with Figure A-1l (Ontario), where there are 9 overlaps at 


minimum and 47 at maximum (opposite the radar). 


A sweep-by-sweep analysis of the Long Beach radar was 
used to determine the overlap environment for aircraft #399. 
The sweep analysis was used to determine how an actual environ- 
ment may interact in producing overlapping replies in space as 
the interrogator scans in azimuth. Figure A-13 summarizes the 
output of the program as it scanned target #399 by 0.225-degree 
sweeps (simulating a four-second rotation rate and 400 pulse/ 
second repetition frequency). The program determines which 
aircraft cause garble to target #399 and which cause potential 
garble. A maximum of eight aircraft replied together, and 
three had overlaps causing serious garble. It should be noted 
that if the system were capable of reading through three over- 
lapped replies, it would have centermarked target #399 correct- 
ly. Although it is probable that passive systems will degarble 
the three overlaps in the Long Beach situation, this informa- 
tion may be unusable, as the radar is too close to the target 
(rapidly changing geometry) and BCAS for smooth tracking. 


Additional plots of the garble contours of the other radars 
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are presented in Attachment 1 in order of decreasing target 


density. 


A.5 STATIC SIMULATION OF LA BASIN, HALF DENSITY 


An analysis similar to that in Section A.3 was con- ‘ 
ducted using the half-density traffic model of 382 targets. 

Figure A-14 shows that the garble contours have been roughly 

halved. Figure A-15 compares the active and passive half- 

density garble levels for various azimuth angles around air- 

craft #399. The ratio of active to passive garble has 

remained at the same 4-to-1 as in the full-density model. The 
half-density Long Beach garble contour is given in Figure A-16, 
showing a significantly improved garble environment. But now 
a tracker that uses only information "in the clear" (ungarbled) 
aS being positively identified would cause cause an error of 


0.5 degree in centermarking the azimuth to #399. Figure A-17 


summarizes program output on a sweep-by-sweep basis. 


A.6 STATIC SIMULATION OF LA BASIN, QUARTER DENSITY 

Again, a complete analysis as in Section A.3 was con- 
ducted for the quarter-density traffic model with 177 targets. 
Figures A-18 through A-21 show these results. It is signifi- 
cant to point out that once again the garble environment has 
improved for the Long Beach radar, yet the centermark error 


has worsened, further degrading its usefulness. A one-degree 
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error in azimuth now exists. This can cause threat range 


errors of up to 0.5 nm and bearing errors of 1.1 degrees. 


A comparison of the three levels of passive garble for 
the three model densities is shown in Figure A-22, indicating 


how the traffic varies with azimuth from aircraft #399. 


A.7 STATIC GARBLE ANALYSIS SUMMARY 

A compilation of the garble analysis of each radar 
for the three model densities is presented in Figure A-23. 
The numbers indicate the maximum passive garble level for each 
radar, which almost never corresponds to a maximum position 
for another radar unless they are approximately aligned with 
target #399. The active garble for each density is also 


given in the last column. 


A.8 ACTIVE GARBLE DIRECTIONAL ANTENNA IMPROVEMENTS 

A static simulation was developed to investigate the 
active garble improvement by employing a directional active 
antenna system. The simulation included antenna directional 
wedge angles of 16° a5", and 30°, which were rotated 360° 
about the BCAS aircraft test position. The BCAS positions 
were chosen to simulate areas where target #399 would be in 


synchronous passive garble and active interrogations would be 
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necessary to augment tracking data. This meant placing the 
BCAS on the opposite side of target #399 from the selected 
radars. Figure A-24 summarizes the statistics of the direc- 
tional improvement over the omni antenna for the three model 
densities. Figure A-25 shows the specific wedges of the 
d.rectional scan and BCAS position which includes both target 


#399 and the radar in a passive garble geometry. 


A.9 DYNAMIC, THREE-DIMENSIONAL SIMULATION 

In order to verify the validity of some of the assump- 
tions made in the formulation of the static, two-dimensional 
simulations, and to allow detailed investigation of some of 
the more complex garble phenomena, a dynamic, three-dimensional 
simulation was developed in Fortran IV for a CDC Cyber 175 


Computer. The program accepts as input a complete description 


of the interrogator-aircraft environment that it is to repro- 

duce. Interrogator characteristics specified for the simula- 
tion include latitude, longitude, scan (rotation) rate, 

interrogation repetition rate (including a provision for six- 

interval and eight-interval staggered interrogation rates), 

and interrogation beamwidth (assumed to be 4° for purposes of é 
the simulations whose description follows). Although the 

garble simulations concentrated on the use of a single inter- 

rogator at one time, some studies of fruit incorporated 70 


ground sites simultaneously—all within 290 nm of LAX. Data 
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Figure A-25 Active Directional Antenna Wedge Garble Toward 
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regarding interrogator characteristics were obtained from the 
file maintained by the Frequency Assignment Staff in FAA's 
Airways Facilities Service. These characteristics are shown 


in the program output listings included in Attachment 2. 


Aircraft deployment made use of the same snapshot of 
the Los Angeles Basin model used for the static, two-dimen- 
sional simulations. The entire group of 743 aircraft was 


processed by the simulation. In addition to x- and y-coordi- 


ae ae Ra ay a 


nates of position, aircraft altitude, climb rate, and horizon- 
tal velocity were provided to and processed by the program. 
The "dynamic" nature of the program is due to its ability to 
account for aircraft motion (including that of the BCAS air- 
craft) as well as for rotation of the interrogators. Turn 
rates specified in the Los Angeles Basin model are ignored; 
instead, the aircraft are moved in straight lines according to 
the horizontal velocity components specified in the model. 

All aircraft are assumed to be equipped with a transponder 
that radiates at a particular power level (nominally 500 watts 
at the antenna end of the transmission line) and with a 
specified reply probability. To minimize the execution time 
of the program, a parameter is required to specify the interval 
at which aircraft positions are updated. BCAS aircraft 
position and velocity in three dimensions are specified ina 


similar manner. 
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The dynamic simulation uses a basic "frame time" of 
2.5 ms within which to compute the various interrogator- 
aircraft interactions. The simulation begins with a random 
selection of interrogator azimuths and interrogation pulse 
start times. Thereafter, the azimuth and interrogation time 
are propagated according to the rate specified in the input 
file. Any interrogation that occurs during a current time 
frame may elicit a reply during either the current frame or 
the subsequent frame. Most of the simulation program logic 
is concerned with the bookkeeping tasks associated with , 


accounting for the interrogations and replies as they occur 


in separate time frames. 


For simulations of garble phenomena, two windows are 
computed. Both involve a protection distance which, although 
parametric within the program, has been assumed to be 20 nm | 
in the computer runs made to date. First, an azimuth window 
is established to ensure proximity among BCAS, target aircraft, 
and the main beam of the interrogator eliciting replies from 
the target. The azimuth window is opened with respect to a 


particular interrogator when the boresight of the main beam 


is separated from the line of sight to the BCAS aircraft by 


no more than the angle 6, where é 


A-38 


In this equation, 
d. = the protection distance (20 nm) 
= the projection of the vector from 
the interrogator of interest to the 
BCAS aircraft in the horizontal 
plane 


Y = interrogator beam width 


Having established an azimuth window, the program computes a 
time window within which replies from a target aircraft must 
fall if the aircraft is within the protection distance. The 
time window begins three microseconds (to account for trans- 
ponder delay) after the interrogation from the ground site of 
interest arrives at the BCAS aircraft, and ends approximately 
250 microseconds later, when the interrogation would have had 
time to travel past BCAS out to the extreme of the protection 
range and back again to BCAS—a distance of some 40 nm. An 
aircraft located at this extreme would have taken the longest 
time to be detected by BCAS, yet still be within the protected 


range. 


Once the two windows are established, the program 
begins an iterative process of reply generation for all air- 
craft/interrogator interactions during the frame. If no 
interrogation occurs during the current frame, then no 


further action is required with regard to that particular 


interrogator during the frame. Otherwise, the interrogation 


time and interrogator azimuth are computed. Next, a test is 


made to see if each aircraft is located within the main beam 


of the interrogator of interest. If it is not, processing 


continues with the next aircraft; if the aircraft is within 
the main beam, another test is performed to ensure that the 
aircraft is visible above the radar horizon. Although the 
aircraft movement geometry is executed for a flat earth, the 
equivalent radar horizon calculation is made assuming a 
spherical earth with four-thirds its actual radius. Finally, 
a random number generator is employed in conjunction with the 
reply probability parameter supplied as input to the program 
in order to determine if a reply will be generated. Once the 
reply probability test is passed, a reply is generated, and 
its arrival time and power at BCAS stored in an array for 


subsequent processing. 


This iterative process continues until all replies 
have been generated (in no particular time order) and stored 
in a reply array. At the end of the frame, the reply array 
is ordered chronologically and made available for printing 
or subsequent processing. The incorporation of reply pulse 
arrival time and power, as well as the inclusion of aircraft 
and interrogator numbers associated with each reply, makes 
further detailed analysis of garble and reply processing a 
straightforward task. Aircraft identification, for example, 


will allow altitude information pulses to be generated when 


er Te ee ee 


specific degarblers and reply processors are ready to be 
evaluated. Even the effects of pulse stretching and phase 
cancellation can be analyzed in detail with such a model. 


These phenomena are the subjects of continuing investigations. 


A.10 RESULTS OF THE DYNAMIC SIMULATION 

The dynamic simulation has been used to validate 
some of the data obtained in the two-dimensional static 
simulations. By placing the BCAS aircraft at a location 10 nm 
south of aircraft #399, for example, the peak number of 
garbled replies was found to correlate with the static model 
to within 13%. Within the windows for the Norton AFB radar, 
39 garbled replies were found in the three-dimensional simula- 
tion versus 45 in the two-dimensional model. One reason for 
the smaller number in the three-dimensional simulation is the 
fact that several aircraft were too low to be interrogated by 
the Norton radar, although they were within the protected 
range of the BCAS. The two models can thus be reconciled 


despite only minor differences. 


Fruit counts obtained for the entire environment of 
70 interrogators and 743 aircraft also correlate well with 
both measurements and predictions of interrogations and 


replies reported previously. Fruit averaged 62,500 replies 


per second and peaked at 84,000 per second over the period 


| 


measured. The lowest rate detected during the period was 


44,000 replies per second. 


Perhaps the most important results obtained to date 
with the three-dimensional simulation will allow important 
decisions to be made regarding radar selection and mode 
selection logic. These results indicate the average and 
maximum number of garbled replies detected on a complete scan 
of the nine radars of primary interest within the basin. The 
BCAS aircraft was placed at 10-nm intervals along radials 
extending north, east, south, and west from LAX. At each 
position of the BCAS aircraft, the number of garbled replies 
was determined for each sweep of the radar (each time window) 
within its azimuth window. The extreme values and the average 
over all sweeps were recorded. One garble is caused by two 
overlapping replies; two garbles by three overlapping replies; 
and so on. A single reply, of course, is ungarbled. No 
reply during a sweep was not included in the average. The 
minimum number of garbles in every case except one was zero; 
that is, in nearly every geometry, including the very worst, 
there was at least one sweep in which some reply arrived in 
the clear. The lone exception had a single garble as its 
minimum. The complete results are presented in Figures 


A-26 through A-29. 


Similar runs were made for the operational 068° 
approach to LAX and for the southern departure at 120°. These 


results are presented in Figures A-30 and A-3l. 
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DISTANCE (NM) FROM LAX ALONG NORTH RADIAL 


RADAR 10 20 30 40 50 60 

LAX 8/9 L.7/21 2.7/11 3.77/12 4.97/14 6.17/19 
Burbank 8/5 APT 3/3 - 7/4 1.0/6 1.5/8 
El Toro 9.2/24 9.8/25 10.7/24 10.6/26 11.6/27 11.8/25 


Long Beach 2.3/9 3-L/1LS., 4:5/17 4.97/17; 65/18 7.6718 


March 11.9/55  9.7/54.. 6.6/43..4.2/13 4.0/11 -3.9/11 
Norton 7.9/41.. 5.6730, 2.39/17 2.57/10. 3.0710 2.9/10 
Ontario 5.6/25 4.1/20 2.3/15 1.8/9 2.3/9 2.4/9 


San Pedro 2.6/13,. 5.2/19.. 7.1/20,. 9.0/21. 10.9/22 12.7/23 


Santa Ana 6.3/21 7.2/24 8.4/23 8.8/23 10.1/23 10.8/21 


Figure A-26 Average/Peak Garble for North Radial from LAX 
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Figure A-27 


DISTANCE (NM) FROM LAX ALONG EAST RADIAL 
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Average/Peak Garble for East Radial from LAX 
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Figure A-28 


DISTANCE (NM) FROM LAX ALONG SOUTH RADIAL 

10 20 30 40 50 60 
1.0/9 1,6/01 2.3742 3.17732 2.2/12 3.8732 
362/15* 6.11/22 $6722 11.8727 11.5/23 LL.T/23 
S20/25. 654/25 4.4/1 2.5710. 1.if? 6/2 
1.5/8 261/78 208/77 33/9 a. 0/9 3.9/9 
$5.6/33 16.12/33 14:1/32 10.47/30 G.if25 3.2/9 
16.23/39 19.3/37 LSst/38 19. 2/37 6.38/37 12.3733 
10.5/24 13.4/27 16.6/27 .16.4/29 15.1/29 13.7/28 
-5/6 safe 33/2 ENLS 8/3 <173 


S.2/15 4.3/16 345/13 -2.8/13 £.0f13 i4/7 


Average/Peak Garble for South Radial from LAX 
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Figure A-29 


DISTANCE (NM) FROM LAX ALONG WEST RADIAL 
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Average/Peak Garble for West Radial from LAX 
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Figure A-30 
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DISTANCE (NM) FROM LAX ALONG 120° RADIAL 


10 20 30 40 50 60 


LAX LPes/9 268/11 2 S972 Tae 2/26 VS. 729 
Buxbank 2.8/9 4.3/16 c 12.2/26.15.0/26 18.3729 
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Figure A-31 Average/Peak Garble for 120° Radial from LAX 
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FIGURE A-32 


LA BASIN GARBLE DATA PLOT 


(SEE DRAWINGS PACKET) 
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APPENDIX A 


ATTACHMENT 1 


The following are passive garble contours for target #399 in 
the LA Basin plotted with respect to each of 9 radars for 


three Basin traffic densities. 
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APPENDIX A 


ATTACHMENT 2 


The following pages contain a source listing and samples 


of output from the dynamic, three-dimensional BCAS 


simulation. 


£2.“ ACAS 


“4 


74/74 or T=2 FIN 4644378 03/35/78 


PROGRAM BCAS CINPUT,OUTPUT,TAPE1,TAPE6,TAPELC) 
C CMMON/ATRCRF /NAC , RHOLJUX (743) ,RHOLIY (743) AL T(743) sRHLUXO(743) 5 —, | 
RHLJYO (743) yAL TOOT (743) »RHLIXC (743) »RHLSY 8 (743) yRNGHCR(743), | 


2L AC (743) y MOSFCW( 743) yXPOW,ALTO (743), ALTSQ(743) 


1 


2 


ao 


cr 
v 


Gets ve 


Ceeeee 


C 


C COMMON/INTER/NINTS gRHOLIX(703,RHOLTY (75), INTIO(70) »ROTRAT(7G), 
PPP(76 48) sINTRVL (70) gLINT(70),AZSTRI (70) ,ISTAG (70) | 
COMMON /BCAS/REOLEX y RHOLBY sAL TB yRH1L3XOyRHLAYID,ALTBOT,RHIBXC,PHIBYC, . | 
ALT80 5RHOUBC (768) »RHOTBC(70) 
COMMON/VE CAG /RHOFCN (768,70) | 
COMMON/TIME/TIME ,TNEXT,FRAMET,ENOTIMyACTIM=, IACUPD yFRAMSC | 
C OMMON/TINEIN/TIMINI (70) »TNXINT(7C),TINSTRK(70) | 


C ONMON/INFUT/COSHBN, REPROB, OPROT yBHOVR2,NTRACK 

C COMMON/WINDOW/IWIN, ITAZWIN(9) sTHNBEG(9,3) y THNEND(9, 3) gp TMNFRT (9) » 

TMXPRIL 

COMMON /CONST/FT2NMyREARTH, THOG3EyPIy HALFP1yTWOPIyOTRyRTDyCyONEOVC 

C OMMON/SORT/SCRT (2990) »ISORT (2020) ySORTNX (2309) ,ISCRTX (2009) ee z 

C CHMON/REPLY /REPNOW (1095) yNOWIUP (1960) »NOWR=P y NEXREF, REPTIM, ISAVE ] 
OIMENSION ISTROL (9) 


ESTASLISH CONSTANTS ANC INITIAL VALUES OF VARIOUS PARAMETERS. READ 
INPUT FILES AND PARAMETERS. 

TIME=0. 

TNEXT=TIME+FRAMET 

ICONT=4 

NFRAME=1 ’ : : ae 

CALL SETUP 
UPDATE BCAS POSITION ANO OETERMINE AZIMUTH WIKCOWS FOR UP TO 9 
INTERROGATORS OF INTEREST.~ 

CALL WINDOh- Seu % ae P= 
IF BCAS IS NOT WITHIN AN AZIMUTH WINOOW, UPDATE TIME AND TRY AGAIN 


ee ae ad 


Cos Ce 0) C¢ 


QO 


eee ne 
v 


fe er ey 


rn 
C% FPtF 
Cte ee 
CeeEre 


Cc 


IF (IWIN. EC.0) GO TO 2 


__ OTHERWISE CHECK FOR THE REGINNING OF A SEQUENCE OF CONSECUTIVE 
FRAMES. IF IT IS NOT THE 3FGINNING, JUMP DOWN ANDO DETERMINE TIME 
WINDOKS FOR RELEVANT INTERROGATORS. IF IT IS THE BEGINNING, 
INITIALIZE REPLY ARRAYS, BACK UP TIME 3Y ONE FRAME, AND DETERMINE 


INITIAL INTERROGATION TIMES. 


IF ({ICONT»EQ.1)? GO TO | 

NOWREP=0 

NEKREP=0 

CALL MEMSET (#165S02T,299C6) 

CALL MEMSET (#1eySORTNX,20)0) | 
NFRAME=1 | 


ON INITIAL FRAME, 30 NOT SACK UP TO PREVICUS FRAME. 
IF (TIME-EC.G5.) GO TO 13 
TNEXT=TIME 
TIME=TIMESFRAPMET 
DETERMINE INITIAL INTERROGATION TINES, INCLUOING STAGGERED PRF. 
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Mt 


BCAS 


rm oD 
GI 


£2 


11 


_ISTRDL CM) =0 


wo 


eens 


S43) Ot 


af 


T4/74 OPT=2 FIN 4e6¢4326 C3/30/78 1% 


CO 7 I=1,NINTS aes 
IF (LSTAG(I) .E%.4) GO TO 6 
TOTPRP=6 

IS=ISTAG(I) 

OO @ K=t,T8 
TOTPRE=TOTPRP #PRE (I,K) 
CONNINUG. as 
TTEMP=TIME FTO TPRP=AHOD(TIME-TINSTR(I),TOTPR) 
IF (TIME*TINSTROI) LT oO0) TTEMP=TINSTR(1) +TIME 
00 9 K=1,IS 
TTEST=TTEMF=FRP(T,ISTAG(I) -K +1) 
IF (TTEST»LTeTIME) GO TO 11 
TTEMP=TTEST _ 

CONTINUE 

PRINT 12 : 
FORMAT (* INITIAL STAGGER TIME ERROR*) 
sToP 

TIMINT (I) =TTEMP 
LINTRVL (I) =ISTAG(T) “K+4 
GO TO 7 

TIMINT (L) =TIMEFPRE (1,1) -AMOO(TIME=TINSTR(I), PRP (1) 1)) 
IF (TIMESTINSTQICI) oh Todo) TIMINT(TYSTINSTR(I) 

C GNTINUE 

O00 29 M=1,NTRACK 


CONTINUE 
* DETERMINE TIME WINDOWS FOR RELEVANT INTERRCGATORS 


O00 15 M=1,\TRACK 

IF (TAZWIN(M).EQ.5) GO TO 15 

IF CISTAG(M).2E9.1) GO TO 22 

ISL=INTRVL (MM) 44 

IF (I1S1eGTeISTAG(M)) ISi=1 

Ts2=Isit1 

IF (I1S2eGTeISTAG(M)) IS2=1 

IF (ISTROL(M) .EQ.0) GO TO 24 +o 23 oe. 
IF (ITSTROL(M) -FEQ61) GO TO 26 
ISTROL (M) =6 

TWNBEG (M,1)=TWNBEG (My2) 

TWNENOD (My1)=TWNEND (M2) 

TWNBEG (HM, 2.=TWNSEG (4,3) 

TWNENG (My 2)=TWNEND (43) 

THNBEG (My Sd=THNBEG (M2) + PRP(M,IS1) 
TVINEND (My 3)=TWN3EG (My 3) + TMXERT 

IF (CTIMINT CM) eh Te TNEXTeANO Oo THNENOD (My 1) GT TNEXTY TSTROLCNIH4 
IF CCTIMINT(M)*PRPEOM,ISL)) oh Te TNEXT*eANDe TWNEND OM, 2) eGToT NEXT) 
LISTROL (M) =2 

GO TO 14 

ISTROL (CM) =6 

ITF (TIMINT (MD) LT oe TNEXT CAND e THNEND (My 1) eGTeTNEXTY ISTROL (MI) =H4 
IF CCTIMINTOM)FPRP(M,ISL)) eL Te TNEXT* ANDe THNEND (M92) GT eT NEXT) 
LISTROL (M) =2 

GO TO 14 

TWNOEG (My L)=TIMINT CM) tTHNPRT (4) 

THNENO (My 1) =TWNIEG OM, 1) (TMXERT 


A=75 


neers 


SuUeo mei i £6 Pobie Tews wevry wveverry * 


THN3EG (My 2)=THNSEG (My 4) #PRP(M,IS1) 
TWNEND (My 2)=TWNOEG (My 2) #TMXPRT 
TWNBEG (My 3)=TWNSEG (My 2) #PRP (My IS2) 
THNENO (My 3)=ThNSEG (My 3) *TMXPRT 
IF (TIMINT (M) eL Te TNEXTeANDe THNEND(M,1) oGTeTNEXT) ISTROL(M) =4 
IF COTIMINT(M)+#PRP (My ISL) oh Te TNEX Ts AND o THNENO (My 2) oGT eT NEXT) . 
LISTROL (mM) =2 
POE Le eee ee 
22 IF (ISTROL(M) «EQ.6) GO TO 2) 
IF (ISTROL(M) E961) GO TO 27 
ISTROL (M) =C 
TWNBEG (My 1)=TWNREG(M, 2) 
TWNEND(M,1)=TWNEND (My 2) 
_ THNBEG (My 2)=THNBEG(M, 3) 
T WNEND (My 2)=TWNEND (My 3) 
THN3EG (M, 3)=TKNSEG(M, 2) #PRP (My 41) 
TWNEND (My 3)=ThNBEG (My 3) #TMXPRT 
IF (TIMINT (M) oLT.TNEXT eAND«s TWNENO (My, 41) 6GTeTNEXT) ISTROL(M)D=1 
IF CCTIMINTOM) #PRP (My 4)) oLT. TNEXT AND. THNEND (My 2) eGToTNEXT) 
LISTROL(M) =2 
CO TO 14 
27 ISTROL (M) =O 
IF (TIMINT (M) oLToTNEXTsANDeTWNENO (My 1) eGTeTNEXT) ISTROL(M) =1 
IF (CTIMINT(M) #PRP(M51)) oL Ts TNEXT «ANDe THNEND (My 2) eGTeTNEXT) 
LISTROL (Mm) =2 
GO TO 14 re 
26 TWNBEG (My 12=TIMINT (M) #TMNPRT (M) 
TWNEND (My 1)=TRNGEG (My 1) #TMXPRT 
THNBEG (My 2)=TWNBEG (M51) #PRP (My 4) 
TWNEND (My 2)=THNBEG (My 2) *TMXPPT 
TWNBEG (My 3)=THNSEG (My 2) +PRP (My 1) 
TWNENO (My 3)=TWNBEG(M, 3) tT MXPRT . : ae 
IF (TIMINT (4) eh Te TNEXT eANDe TWNEND (My 1) oGTeTNEXT) ISTROL(M) =i 
IF ((TIMINTOCM)*PRP (My1)) oh Te TNEXTs AND oe THNEND (My 2) oGToTNEXT) 
LISTROL (M)=2 
14 PRINT 114 INTIO(M) yTIMEyTNEXTyTIMINT (M) y (CTHNBEG (My K) y THNEND (MyK) y 
1K=1,3) 
44% FORWAT 6/ F295 F3; SF 11.6 752% SNINY 3B C2Keeeea eee) tis 
15 CONTINUE - 


CALL FRAME 
ICONT=4 
nena ORDER REPLIES CHRONOLOGICALLY SY PACKING PRE=SORTED ARRAY ANC 
reese APPLYING 3UB3LE SORT if S20 ee 


If (NOWREP.EC-9) GO TO 18 
CAlLt CHRONC 


PRINT RESULTS 
‘i yreut 


Owe @EPLY BRAY F2OM NEXT TO CURRENT FRAME CUNPACKED) 


sé 5 


74/74 oP T=2 FIN Ged t4o5n 


bor ourtiy aor 


CALL MEMSET (16,SORTNX, 26635) 


senee UPDATE AIRCRAFT POSITICNS (AND RELATED PARAMETERS) AT APPROPRIATE 
sewer INTERVAL. 


ro TE Ms he 


IF (NFRAME*EGeIACUPD) GO TO 5 
NFRAME=NFRAME #1 

60.16 3 
NFRAME=1 
PRINT 105,RHO18X sRHOLBY, ALT,» (JyRHOL IX (J) pRHOLIY (J) ,AL TCU), 
LRHOJIC(J) »J=4 NAC) : 

FORMAT (1Hi,s*@CAS POSITION IS *3F13.2///* AIRCRAFT DATA (NOe 
£ Y ALT RANGE FROM BCAS) *//( IG 2F8e2 yp FIel gy F702 gl hX Ib y2FBe2y 
IEG Ge Foe sIeX Tap SEOe Osh Sw A ghT eA). enc yng sien tea 
CALL MOVEAC 

GO TO 3 

ICONT=0 


ann ne UPOATE TIME ANDO CHECK FOR END OF SIMULATION. 


WDOOON 


TIME=TNEXT 

IF (TIME-GTsENOTIM) GO TO 10 
TNEXT=TNEXT+FRAMET 

GO TO 1 f 
ERKOFILE 10 


REWIND 10 | 
STOP 
ENO 


RELOCATION j Wis 
TIME REAL AIRCRF 
BCAS REAL BCAS 
BCAS REAL AIRCRF 
ALRCRF REAL AIRCRF 
INTER [ REAL ae -INPUT 
CONST REAL INPUT 
INPUT REAL CCNST 
TIME e REAL TIME 
TIME REAL CCNST 
CONST INTEGER 
TIME INTEGER WINDOW 
INTEGER INTER 
INTER INTEGER 
REPLY INTEGER SCRT 
SORT INTEGER IATER 
INTGeP INTEGER 
INTEGER INTEGER WINDOW 
INT*=GER INTEGER 


20UTINE WINDOW 74/74 OP T=2 FIN 46644338 GI/357% 


SUBROUTINE WINDOW 

C OMMON/BCAS/RHOL BX y HOLBY y ALTB y RH1 BX Dy RH1BYD yALTEDT yRH18XUymHIBYL, 
1ALT89 yRHOJEC (7638) »RHOIBC (76) 

C OMMON/WINCOW/IWINy LAZHIN(9) »TWNBEG(9, 3) y THNEND (9, 3) > TNNPRT (9) 4 
1 TMXPRT 

C OMMON/INT ER/NINTS yRHO1LIX (72) yPRHOLTY (70) yp INTIOC7G) yROTRAT(7G) 5 
LPRP (76 98) y INTRVL (79) yLINT(73) 5 AZSTRI (70) y ISTAG (76) 
__ COMMON/TIME/ TIME, TNEXT yFRAMETy ENDTIM yACTIMEy LACUPD yFRAMSC | 

COMMON /INFUT /COS HBWy REPROB, CFROT y BHOVR2 yNTRACK 

C OMMON/CONST/FT2NM,REARTH, THOS 3£,PI,HALFPI yTWOPI yDTRyRTD 5C yONEOVE 

C OMMON/AZ/AZIAT 

IWIN=0 


a 
v 


Cerrs COMPUTE BCAS POSITION 
Cc 


-RHOLBX=RH1 EXO +RH1BXO* TIME 
RHOLBY=RHLEYO #RH1BYO* TIME 
ALTB=ALT904AL TBOT* TIME 
00 10 I=1,\TRACK 

___ TAZWIN(T) =6 

( 

Crenee COMPUTE INTERROGATOR AZIMUTH 

c 
AZINT=AZSTRT CI) *ROTRATCL) *TIME 
AZX=SINCAZINT) 

AZY=COS(AZINT) 
RHOIBX="RHC1IX(T) #RHOLBX 
RHOLBY==RHC1TY(I) *RHOLBY 
RHBXY2=RHOIS8X*RHOI 8X +RHOLBY *RHOIBY 
RHOISC (L) =SQRT(RHSXY 2+AL TB*ALTA) 
RHBCXY=SQRT(RESXY2) 
_TMNPRT (I) =RHCIBC(ID/C + 3.E*06 
CPRTAN=COS (ATAN (DPROT/RHBCXY) #8WOVR2) 
C OSANG=(AZX*RHOI BX+AZY*RHOIBY) /RHBCXY 
IF (COSANGeLT.CPRTAN) GO TO 19 
IAZWIN(I)=4 
IWIN=4 

13. CONTINUE _ 
RETURN 
END 


FCLIC <xEFERENCE MOP (R=1) 


SN TYPE RELCCATION 
evs QEAL TIME 2 ALTB REAL BCAS 
eo REAL BCAS 19 ALTBQ REAL BCAS 
tT REAL AZ 1725 AZSTRT REAL ARRAY INTER 
REAL 63 AZY REAL 
as REAL INPUT tc 6C¢ REAL CCNST | 
his REAL 9 COSHBW REAL INPUT 


7TH/74 O° T=2 FIN 4e64+4338 Q3/30/78 


SUBROUTINE FRAME 
C OMMON/ALQC2F/NAC y RHOLUX (743) y PHOLIY (743) y ALT (743) yRHLIXD (743) y 
LRHLJSYU (743)  ALTOOT (743) y RH1IXE (743) RHLIYG (743) pRNGHCRITZ43) y 

2L AC (7&3) yMCBFCW( 743) sXPOW,ALTO (743) y ALTSQU743) 

COMMON /BCAS/RFOLAX yRHOLBY yALTByRH1 9X Dy RH1 BYD,ALTBOT,RH1IBXA2RHIBYC, 
LALTB% yRHOJEC (768) »yRHOIBC (70) 
CGMMON/INTER/NINT Sy RHOLTX (72) » RHOLTY (70) py INTIO(70) sROTRAT(7L), 
LPRP (7G 98) yINTRVL (79) gLINT (73) sAZSTRT (705 ISTAG (70) 
C OMMON/VECMAG/RHOFCN (768470) 

C CMMON/TIME/T IME y TNEXTyFRAMET, ENDT IM ,ACTIMEy IACUPC yFRAMSC 
COMMON/TIMEIN/TIMINT (70) yTNXINT (70), TINSTR(70) 

C OMMON/INPLT /COSHAW,y REPROBy OPROT BHO VR2yNTRACK 

C OMMON/CONST /FT2NM, REARTH, THOG3E ,PI,HALFPT,TWOPI yDTRyRTD yCyONEOVCE 

C OMMON/SORT/SORT (2000) ,ISORT (2090) ySORTNX(2200),ISORTX(2060) 
COMMON /RE PLY/REPNOW (1009) yNCHIJP( 1030) yNOWREPNEXREF »REPTIMy ISAVE 


Cena START OF MAIN INTERROGATION LOOP : 
, 50 15 I=1,NINTS 

cress IF THERE IS NO INTERROGATION THIS FRAME, GO TC NEXT INTERROGATION} 
‘ IF (TIMINT(I).GEsTNEXT) GO TO 15 

alia UPOATE NEXT INTERROGATICN TIME, ALLOWING FCR STAGGERED PRF. 


IF CISTAG(I) -EN.1)? GO TO 2 
INTRVL CID =INTRVLCTI) +1 

IF CINTRVL (I) GToISTAG(I)) INTRVL CT) 1 
TNXINT (IV =TIMINT (I) *PRPCI,INTRVL(TI)) 


iaaiaea SET SECCND FRAME INTERRCGATION FLAG TO APPROPRIATE VALUE. 


oOaan 


INTNO2Z=4 
IF (CTNXINT(L) LTeTNEXT) INTNOZ=2 
c 
Crees CALCULATE INTERROGATOR PEAM AZIMUTH 
c i Coxe eM 
AZINT=AZSTRT (1) *FOTRAT(I)*TIME 
A 7X=SINCAZINT) 
AZY=COS(AZINT) 


. 
Cush ee COMPUTE INTERMEDIATE VARIASLES 
c 


TINTP3=TIMINT (1) +36&°06 ‘ og 
TNEXP3=TNXINT (1) +326 *U6 


C iI 
Catal os START OF MAIN AIRCRAFT LOOP 
c 
00 12 J=1,NAC 
Sith CALCULATE DCT PRODUCT OF INTERROGATUOR AZIMUTH UNIT VECTOR ANC RANGE 
Ceerey VECTOR FROM INTERROGATCR TO AIRCRAFT. 


lol 
u 


CTHETA=(AZX* (RHOLIX (J) SRHOLIX(T)) tAZY*# (CRHOLIY (JP eRECLTY(I)))* 
LRHOFCN (J,T) 


FRAME 74/74 OPT=2 — FIN 4.644328 93/30/78 1 


cease IF AIRCRAFT IS NOT WITHIN INTERROGATOR MAIN BEAM, GC TO NEXT AIRCRAFT 
; IF (CTHETAsLTeCOSHBH) GO TO 10 Li aie Teyts A aa lee 
Ganees IF AIRCRAFT IS NOT VISIBLE ABOVE THE HORIZON, GO TO NEXT AIRCRAFT 
: RHOISM=SORT((40/ (RHOFCN( J, I) *RHOFCN(JyI))) HALTSC(J)) 

IF (RHUIUMeGEsRNGHOR(J)) GO TO 10 


C 
Czeeer IF FIRST REPLY IS NOT MISSED OUE TO REPLY PROEABILITY, CHECK FOR A. 
Green SECOND INTERROGATION THIS FRAME. 
c 

IF (RANF (DUM) LE*REPROB) GO TO 4 
Cc ooo wees ee oe - — eee tee: we <ianirenigareanimrermeneneeerin — SS ae 
guaree FIRST REPLY MISSED DUE TO REPLY PROSASILITY. IF NO SECOND 
cust INTERROGATION THIS FRAME, GO TO NEXT ATRCRAFT« 
c 

IF (INTNO2.EQ.1) GO TO 10 
c 
cases "OTHERWISE, CHECK FOR SECOND PEPLY PROBABILITY. IF MISSED, GO TO NEXT 
Cretes AIRCRAFT. 
c 

IF (RANFC(OLM) .GTeREPFROB) GO TO 10 


rn 
uv 


Caress OTHERWISE, COMPUTE ARRIVAL TIME FOR SECOND REFLY AND STORE If 
c#eeee == APPROFRIATE ARRAYs © 
ce 


7) 


REPTIM=TNEXP3 #(RHOIJM#RHOUBC (J) ) *ONZOVC 
IF (REPTIM+GEsTNEXT) GO TO 3 
ISAVE=LINT (I) *MDBPOW (J) 
CALL SRTNOk 
GO TO 10 + 
3 ISAVE=LINT (I) 4MOBPOW (J) 
CALL SRTNEX 
GO TO 10 


Cneees FIRST REPLY IS OKs LOOK FOR A SECOND INTERROGATION THIS FRAME. 
b IF CINTNO2sEGe1) GOTO 8 AGT e TT rae +e eee are eS 
C 
Crewe SECONC INTERROGATION OCCURS LATER THIS FRAME. COMPUTE FIRST REPLY i] 
Ceeean ARRIVAL TIME. 
G | 

REPTIM=TINTP3 #(QHOTUM#RHOJAS (JI) *CNEOVEG _*®@ 

IF (REPTIM.GEsTNEXT) GO TO 7 
| 
! 


r 
¥ 


Gree ee FIRST REPLY ARRIVAL OCCUFS THIS FRAME. STCRE ACCORDINGLY. 

C 
ISAVE=LINT (1) #*MOBPONW (CU) 
CALL SRTNOh 


c 
v 


Cs de ig CHECK FOR SECOND REPLY FROBABILITY. IF MISSEC,y GO TO NEXT AIRCRAFT. 


fal 


IF (RANF COLUM) -GTeREPFOS) GO TO 10 


5 tana OTHERKISE, COMPUTE ARRIVAL TIME FOR SECOND REFLY ANC STORE IN 
Caen as APPROPRIATE ARPAY. 


A-80 


é 
E 
5 
, 


—ae lle 


FRAME 74/74 oPT=2 FIN 42€#+4338 03/30/78 12. 


o 


REPTIM=REPTIM*PRP(I,INTRVL(T)) 
IF (PEPTINGETNEXT} GO TO 5 
ISAVE=LINT (I) #*MORPOW (J) 

CALL SRTNOW 


GO TO 10 
5 ISAVE=LINT(I) *MOBPOW (J) 
CALL SRTNEX BA Pete SOUT NT 2 oe 
GO TO 10 
c 
caseee FIRST OF TWO REFLIES OCCURS IN NEXT FRAMEs’ STORE ACCORDINGLY~ 
2 
7 ISAVE=LINT (I) *MORPOW (J) 
CALL SRTNEX  —__ id wagh Vaan a ern eee ie ee eae 
cereee CHECK FGR SECOND REPLY FROAASILITY. IF MISSEG, GO TO NEXT AIRCRAFT. 
IF (RANF (COLUM) .GTeREPROS) GO TO if 
c 
C+s¥ee = OTHERWISE, COMPUTE ARRIVAL TIME FOR SECOND REPLY ANC STORE IN NEXT. 
gaaeee FRAME ARRAY. 
c 
REPTIM=REPTIM#PRP(I,INTRVL(I)) 
ISAVE=LINT (I) #*MDBPOW(Y) 
CALL SRTNEX 
GOTO 160. a 
c 
Create ONLY CNE INTERROGATION OCCURS THIS FRAME» COMPUTE RESULTING REPLY 
cesses ARRIVAL TIME AND STORE ACCORDINGLY. 
3 REPTIM=TINTP3 +(RHOTUM#RHOURC (J) *ONEOVC 
IF (REPTIMLT.TNEXT? GO TO 9 
ISAVE=LINT(I) #MOBPOW( J) 
CALL SRTNEX 
GO TO 10 
5 ISAVE=LINT (I) #*MOBPOW (J) 


CALL SRTNOk 
12 CONTINUE 
Cc : 
C*¥e¥e UPDATE IATERROGATION TIME 


re 
~ 


IF (INTNO2eEQe1) GO TO 11 
IF (ISTAG(I)EN.1) GO TO 12 
INTRVL CL) =INTRVL OT) +1 
IF CINTRVL(LT) GTeISTAG(I)) INTRVLCID =4 
Y TNXINT (I) =TNXINT (TY *#PRPCI,INTRVLOI)) 
1 TIMINT CI) =TNXINT (T) 


one ae IF REFLY ARRAY EXCEEOS 940 ELEMENTS, PRINT WARNING. 


IF (NOWRES .GT.990) PRINT 101 
oe | FORMAT (* WARNING = REPLY ASRAY (PEPNOW) EXCEEDS 9fc ELENENTS.*) 
c ONTINUE 

RETURN 

ENO 


setup 74/74 apt=2 FIN 4.644338 C3/3i/78 4 


£ 


260 


alt 


_COMMON/CONST ZFT2NM sREARTH, THOG3E PI, HALFOI yTWOPI yDTRyRTD yy ONEOVC 


SUBROUTINE SETUP 

C OMMON/ALIRCREF /NAC 5 RHCLIX C743) gPHCLIY (743) ,ALT (743) pRHLISXO(743), : . 
ALRHLJV0 (743), AL TOOT (743), RH1 JX 0742), RHLIV I (742) sRNGHOR(743) , 

2LAC (743) sMCBPOW( 743) »XPOW,ALTO (743), ALTSQ(743) 
COMMON/SCAS/RHOLBX »RHOL3Y »ALTPyRH1 3X0 yRHLBYD,ALTBDT yRHIBXOyRHIBYO, 
LALTB9 »RHOJEC (768) yRHOIBC (70) 

C OMMON/INTER/AINTS,RHOLIX (72) ,RHOLTY (70), INTIO(79) yROTRAT(70), 
__1PRP(70 8) gINTRVL (79) yLINT(79) ,AZSTRI (70), ITSTAG(70) 

COMMON /WINCOW/IWIN, TAZWIN(9) yTWNREG(9, 3) THNEND(943) yTMNPRT(9), 
LTHXPRT 

COMMON/TIME/TIME ,TNEXTyFRAMET,ENOTIM,ACTIME, IACUPC yFRAMSC 

C OMMON/TIMEIN/TIMINT (70) »TNXINT(70),TINSTR(70) — 
COMMON/INFUT/COSH BW, REPROB , OPROT , BWOVR2 yNTROCK 


C OMMON/SORT/ SORT (2396) pISORT (2000) ySORTNX (2900) yISCRTX(2000) 

COMMON/REPLY/REPNOW( 1000) pNCWIJ° (1000) ,NIWREP, NEXREF,REPTIM, ISAVE 

DIMENSION LAT(7543)yLONG (75,43) 5CAR(3) ,AQUT (3) 

NAC=743 

FT2NM=10/6076 012 
REARTH=34376747_ 

TWO43SE=2.* (40/36) *REARTH 

PI=3.1415926535898 

HALFPI=6..5*PI 

THOPI=2.*PI 

OTR=PI/18%. 

RTD=1800/P1. 

C=161875. 

ONEOVC=1./C 

READ 1,FRAMET ,ENOTIM,NTRACK, IACUPD,REPRO3, SMWOTH,OCISCINy XPOW , 
1DPROT » RH1 8X3 yRHLBYO, ALT30,RH13X0yPH18Y0,ALTIOF 

FORMAT (2F416-6,215,5F10.C/6F 10 69) 

PRINT 20d, TIME,FRAMET yENDTIMyRH1BX0, RHLBYI »ALTBCyRELEXO,RH1LBYD, 
1 ALTBDT yNACyIACUPD, XPOW,REPROB,DISCIN ,BMHOTH 

PRINT 201,NTRACK yDPROT 

| 


FORMAT (152,*@CAS SIMULATION INPUT PARAMETERS*#/////6H *****,*TIME* 
16H ¥*#***//T5S7,*START TIME =*File7/TS7y*FRAME TIME =*F1ile7/T59, 
2*ENO TIME =*Fi1e7///6H *****,*GEOMETRY* 0H *¥F**// 

3* THE CARTESIAN FRAME USED THROUGHOUT THIS FLAT*EARTH SIMULATION I 
4S CENTERED AT INTERROGATOF NOe 1 WITH X EAST, Y NORTHy ANO Z UPe*/ 
5* A CALCULATICN IS PERFORMED TO DISCOUNT AIRCRAFT THAT WOULO BE Ov 
6ER THE RADAR KORIZON IN A SFHERICALSEARTH MIDELe*///6H *****¥*BCAS* 
76H #¥**#//* ACAS AIRCRAFT PCSITION CNMyNM,FT) AND VELOCITY (KT,KT, 
BFPM) AT TIME=U ARE*/* X =*F7o2y*%y Y =*F7e29*5 Z =*FBaly5*, xO 
JOT =*Fely*,- YOOT =*FGel 5*,y ZOOT =F¥F7olyto*///EH *¥***,* TARGE 
1T AIRSRAFT*,6H *¥****//I4,*% TARGET AIRCRAFT POSITIONS WILL BE UPDAT 
2EO EVERY*1I3,* FRAMES. (POSITIONS ARE LISTED ON SUBSEQUENT PAGES. 
2)*/* ALL AIRCRAFT ARE ASSUMSO TO BE EQUIPPED WITH TRANSPONDERS THA 
4T RADIATE *F4.9,* WATTS OF FCWER WITH NO ATTENUATICN OTHEF THAN A* 
S/* RANGE=DEPENDENT PROPAGATICN LOSSe THE TRANSFONCERS HAVE A PEPL 
6Y PROBABILITY OF *F4a2y* et ///EH FFFE*SFINTERZROGATCRS* 6H FE*E SS 
7* UP TO 70 INTERROGATORS WITHIN *yF4eu,y* NM OF INTERROGATOR NOe 1 
8ARE INCLUOED IN THE SIMULATICNse THEIR INTERRUGATICN BEAMWIOTHS*/ 
B* ARE *F3.1,* DEGREES, AND THE REST OF THEIR CHARACTERISTICS ARE C 
LESCRISED ON TEE NEXT PAGE.*///) 

FORMAT ¢ 6H #¥* #4 FREPLY WINCCWS* »6H F444 // 
2% THE FIRST *11,* INTERROGATORS ARE TRACKED TO OETERMINE AZIMUTH A 
3NOQ TIME WINOQWS WITHIN WHICH AIRCRAFT REPLIES ARE EXAMINED. */ 


A-82 


+ 4/74 OPT=2 FIN 4.644338 03/30/78 134 


SETue 


&* THE CALCULATION OF THESE WINDOWS ASSUMES A PRCTECTEOQ RANGE OF UF 
5 TO *F4e0,* NM FROM THE BCAS*EQUIFPED AIRCRAFT e*) : 
THXPRT=2-*CFRCT/C 
IF (NTRACKeLEe9)? GO TO 9 


PRINT 109 
169 FCRMAT (* ATRACK EXCEEOS 9*) 
STOP ; 
S) CONTINUE = ee Sd Bap ee ee ee eS ee eee ee 


RH1iBXD=RHLEXO/36LO. 
RHI BYO=RHLEYD/3600. 
RHOLAX=RHLEX9FRH1L8X0* TIME 
RHO18Y=RHLEYO+RH1IBYO* TIME 


crts** CONVERT BCAS VELOCITY INPUTS FROM KTS TO NM/SEC (ALTITUDE FRCH FFM 
Er eee+ TO NM/SECe 


ALTBOT=ALTEDOT*FT2NM/60. 
ALTBO0=ALTSUtFT2NM 
ALTB=ALTBOtALTBOT*TIME 
CALL PAGE Es Piet a 
PRINT 16 
16 FORMAT (T3y*I*yT7y*#ID* yT1iby *LAT* 9 T25 5 *LONG*y T 34g *SCAN* 5743 y*AZ(0)* 
151T53,*BANG (9) *,T72,* PULSE REPETITION RATE*,T1I3,*REOLIX*,1T113, 
2,*RHOLIY* »T124,*RANGE*/) 
I=1 
£7 READ 2,INTIOCIY, (LAT CI, K) pK=45 3) 5 (LONG(I5L) yL=193) yROTRAT(ID » 
1 (PRP (I 4M) yM=155) 
2 FORMAT (110,1X_313,1X53I13,F10.0,8F5.0) 
IF CINTIO(I)-EQ-6) GO TO 3 
XLATL= (33 2455 0/602+570/3600.)*0TR 
XLONGL=ATWOFI$ (118 6+24./600+230/36900)*DTR 
CALL LL2CAR (XLAT1,XLONG1,C.e yCAR) 
ANG=HALFPI+*XLCNG1 
CALL ROTZ (CAR,ANG,AOUT) 
ANG=HALFPI*XLATL 
CALL ROTX (AOUT,ANG,CAR) 
R1X=CAR(1L) 
RAY=CAR(2) ah : - Shee le ee ae ee a 
XLAT =(FLOAT(LAT (I 91)) FLOAT ILAT (I 92)) S60 0+FL CAT (LAT(I93))/36006 
1) *DTR 
XLONG | =ThOPI@( FLOAT (LONG (I 9199 #FLOAT(LONG(Iy2))/EC ot 
LFLOAT (CLCNG(I53))/36C0Ce)*OTR 
CALL LL2CAR (XLAT,XLONG,00,CAR) 
ANG=HALFPI#XLCNGL 
CALL ROTZ(CAR»,ANG,AOUT) 
ANGSHALFPI=-XLATA 
CALL RCTX(AOUT,ANG,CAR) 
19 QHOLIX (I) =CAR (1) ©R1X 
RHOLIY (1) =CAR (2)°R1Y 
RNGINT=SART(RHOLIX(T) **2*RHOLTY (I) #*#2) 
IF (RNGINT.GE-DISCIN) GO TO 17 
INTRVL (1) =1 
TINSTR (1) =RANF (DUM) *FRAMET 
TIMINT (I) =TINSTR(I) 
LINTCI)=INTIO(I)* (2** 30) 
AZSTRT (1) =FANF(OUM)*360. 
PRINT 18, I,INTIO(I) y (LAT (I yK) 5 K=453) » (LONG (LT yL) yL21y3) yROTRAT(I) y 


A-83 


Fes a ee ee ee 


A 


13 


606 TO 11 


13 


103 


+ 


__FRAMSC=19990/FRAMET 


. 
L 


Cree 


€ 


™ 
™m 


~IFERE 


TYPE 
Ad 1” 
ZEAL 


y= SETUP 74/7h OPT=2 FIN 4e6+43728 G3/23/78 15 


LAZSTRT (I) »TINSTR(I), (PRP (19M), M=153) sQRHOLIX(I) pRHOLIY (ID sPNGINT 
FORMAT (1X 9I2e2X%5I3,)2Xy3I351X9 1X9 FT3 5 3X pF holy 3X _2K pF 5019 3X yi X yg F Be E 
LeiXs8 F500 2X9 F601 92X 92Xy F601 y2Xy2X F601) 

ROTRAT (CIP =KOTRAT(IS* THOPI/63. 
AZSTRT (I) =AZSTRATCIO*OTR 

IF (PRP(I,Z) -EQeCe) GO TO 1) 
ISTAG(I)=8 

IF (PRPCI,7) oEQeCe) ISTAG(ID=6 pS ee one oe eee 
ISTAG(D)=1 

IS=ISTAG(I) 

00 13 m=1,1S 

PRP (I, MP =1Le/PRP( IM) 

CONTINUE __ 
T=I+1 

IF (IeLE.76) GO TO 17 

PRINT 163 

FORMAT (//€H *¥*¥***,# WARNING = INTERROGATOR LIMIT CF 70 HAS BEEN R 
LEACHED -*//) 

@INTS2I-1 ss 
IF (NINTS-CTeU) GO TO 5 
PRINT & : 
FORMAT (*4180 INTERROGATORE SPECIFIEI.*) 
STOP 

ACTIME=FLOAT (IACUPD) *FRAMET 


BWOVRZ=0TR*C e S*BMWOTH 
COSHBW=COS (8WCVR2) 


* TEMPCRARY REOUCTION TC 1 INTERR0GATOK (FIRST) 


NINTS=1 

CALL BASIN 

X POW=XFOW* (C/ (40 *1990 E+ 06*FI) I**2 

DO 22 J=1,NAC 

RHOJBC (J) =SORT (CRHOLEX*RHOLI XC) ) #*#2 4+ (RHOLBY=RHOLJY (J) ) ¥¥2+ 
1 (CALTB<ALT (J) ) #*2) 

LAC(J) =3* C2" F209) Fs - Be an ee 
MOBPOW (J) =A3S (3009641000 o*ALOGI0O (CX POWs (RHOUIC (J) **2))) 
IFCIABS(MDEPOW(S) ) GT o(2**2I91)) MDSPOH (J) =ISIGN(2**20°1,MDBPOW (J)? 
MNAPOW (J) =MOBPOW (SD +LAC(U) F 

C CNTINUE 

RETURN 

END 


NCE MAP (R=1) 


FEL CCATION 
TIME 2fic | ALT REAL ARRAY AIRCRE 
BCAS 5 ALT830T REAL BCAS 
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SASIN THIS 7TH OPT=2 FIN 42644328 03/30/78 12. 


SUBROUTINE BASIN 
‘ C OMMON/AIRCRE/NAC y RHGLUX (743) yRHOLIY (743) gALT(743) pRHLIXD(743) 4 
LRHLJYO (743), ALTOCT (743) y RHI JXO (743), RHLJYC (743) yRNGHOR(743) y 
2L AC (743) yMCBPCH( 743) »XPOW ALTO (743), ALTSQ (743) 
COMMON/INTER/NINTS yRHOLIX (79) »RHOLIV (70) y INTIOC70) yROTRAT(7U), 
1PRP (75598) yINTRVL (70) »LINT(79),AZSTRT (70) »ISTAG (70) 
. C OMMGN/VECMAG/RHOFCN( 768,470) 
C CMMON/TIME/TIMEgTNEXT»FRAMETyENOTIMyACTIME,TACUPGyFRAMSC 
C OMMON/CONST/ FT2NMyREARTH, ThC4 SE pF Ly HALFPI y» THOPI yCIRyRTD Cy ONEOVE 
00 5 J=1,NAC 
READ(1,2) REL UXO(S) ,RHLSYO (J) gALTO (J) gRHLSXOC I) yRHLUYO (J) gALTOOT (CU) 
2 FORMAT (12X54 (1X E146 7)/20E L407 1X)? 
READ (1) 
ALTO(J)=ALTI CJ) * FT2NM é 
ALTOOT (J) =ALTOOT (J) *FT2NM 
ALT(J) =ALTO(J)#AL TOOT (Ud * TIME 
ALTSO(J)=ALTCJPFALT OUD 
RNGHOR (J) =SORT(THOG3E*ALT (J) 
RHOLSX (J) =RHLIXO (J) *RHLIXD (J) * TIME 
RHOLSY (J) =RH 1 JYO (J) #RH1 IYO ( J) * TIME 
3 DO 5 I=1,NINTS 
RHOFCN (Jy T)=4 o/SORT ( (RHOLSX (J) @RHOLIXCLI) *¥2 
1+ (RHO1 JY (J) =RHOLTY (1) **2) 


5 CONTINUE 
RETURN 
ENO 
REFERENCE MAP (R=4) sa ee eee | 
| 
' TYPE RELCCATION 
oJ ____trHe OTE OT ER RRA eter 
REAL ARRAY ALPCRE 21326 ALTSQ REAL ARRAY “AIRCRF 
REAL ARRAY AIRC2F 14725 AZSTRT REAL ARRAY INTER 
25AL CONST 6 OTR REAL CONST 
REAL TIME 2 FRAMET REAL TIME 
2 SAL TIME 0 FT2NM REAL CONST 
2EAL : CONST _ £67 INTEGER 
INTEGER TIME 215 INTIO INTEGER ARRAY INTER © 
IMITEGER AKRAY INTER 2633 ISTAG INTEGER ARRAY INTER 
INTEGER 45040 LAC INTEGER ARRAY AIRCRF 
INTEGER ARRAY INTER 16497 MNBPOW INTEGER ARRAY AIRORE 
iNTEGER ALFCRF 0 NINTS INTEGER INTER 
REBL CONST Jy PL REAL CCNST 
RESL ARRAY INTER 1 EARTH REAL CCNST 
REAL ARRAY VECMAG 1 RHOLTX REAL ARRAY INTER 
2 EAL ARRAY INTER 4  RHOLUX REAL ARRAY AIRCRF 
QE AL ARRAY ALPCRE 4266 RHLIXD REAL ARRAY AIRCRF 
REAL ARRAY ALFC2F 5635 RH1LJYO REAL AFRAY AIRCRF 
ZEAL ARRAY ALRCRF 42471 RNGHOR REAL AFRAY AIRCRE 
2 EBL ARRAY INTER 7 RTO REAL CCNST 


A-85 


Tus MOVEAC 74/74 OPT=2 FIN 466¢+4358 63/33/78 


SUBROUTINE MCVEAC 

COMMON /AIRCRF/NAC » RHOLIUX (7423) y PHOL SY (743) ALT (743) yPHIUXO(743) g | 
tRHLJYO (743), AL TOOT (743) »RHLIX6 (742), RHLSV0(743) »RNGHOR(743), 

2L AC (743) yMCBPOW( 743) yXPOW,ALTC (743), ALTSQ(743) 

COMMON/BCAS/RHO1 BX, RHOLSY sALTB RH: 3XO,RHLAYI,ZALTBOT,RHIBX2,RHIBYG, 
LALT BO yRHOJBC (768) »yRHOTBC (70) 
COMMON/INTER/NINTS ,yRHOLTX (799 ,RHOLTY (70), INTIO(70) »ROTRAT(70), 


COMMON/VECVYAG/RHOFCN( 768,79) 
COMMON /CONST /FT2NMy REARTHyTWOS3EyPT,HALFPL,TWOPI yOTRyRTO 9C yONEOVC 
90 & J=1,NAC : 
QHOLIX (J) =RH1LIXG (J) #RHLIXD OS) * TIME 
RQHOLJY (J) =RHLUVO (UP FRHLSYOCUDETIME 
IF (ALTDOT(J) oEQ600) GO TO 3 
ALT (J) SALTO(JDFALTDOT (J) * TIME. 
ALTSO (J) =ALT(J)*FALT(U) 
RNGHOR(J) =SQRT(TWO43E*ALT (J) ) 
IF (CALT(J) LEoe01) ALTOOT(J)==ALTOOT (J) 
3 CONTINUE ae : Ge eet Sth fe < Us >) 
RHOUBC (J) =SORT CC PHOLBX=RHOL IX (J) ) **2 + (RHOL SY =RHCLSY (5) *¥ 24 
1 CALTS*ALT (J)) *#2) 
MOSPOW (J) =ABS (3000641000 -*ALOG19 (XPOW/ (RHOS3C (J) **Z) 0) 
IF CIA3S (MCEPOW(S) ) GT. (2**20"1)) MCS°OW (J) =ISIGN (2**20"4)MOSPOW (J)? 
MDBPOW (J) =MOBFOW (J) LAC (J) 
& GONTINUE 
00 5 I=1,NINTS 
O00 5 J=1,NAC 
RHOFCN (J 1)=16/SORT ( (RHOLUX (J) -2HOLTX (1) ) ¥*2 
1+ (RHOLJY (J) -RHOL ITY (19) ** 2) 
5 C CONTINUE 
RETURN 
ENO 


| 
{ 
' 
t 
j 
4 
{ 
i] 
H 
| 
| 
| 
{ 
1m owaltal pense niamnahins Brion 5 asnamen reste sene mere weenie 


410 XE FERENCE MAP (R=1) 


3M TYPE RELCCATION ; 

2EAL TIME e7ir WALT REAL ARRAY AIRCRE 
SEAL BCAS 5 ALTBDOT REAL BCAS 
ZEAL PCAS 720% ALTOOT REAL APRAY AIRCRE 
REAL ARRAY AIRCRF 17757 ~ALTI REAL ARRAY AIRCRE 
REAL ARRAY INTER qa OC REAL CCNST 
QEAL CONST 2 ENOTIM REAL TIME 
SEAL TIME 6 FRAMSC REAL TIME 
RELL CONST & HALFPI REAL CCNST 
INTEGER 5 TfACUPO INTEGER TIME 
INTEGER ARRAY INTER 1541 INTRVL INTEGER ARRAY IATER 
INTEGER ARRAY INTER wire J INTEGER 

INTEGER ARRAY ALRCOREF 4617 LINT INTEGER ARE AY INTER 


INTEGER ARRAY AIRCRF 9 NAC INTEGER AIRCRF 


ALAA LR ACA A 


OQUTINE SRTNOW TH/ 74 O°T=2 FIN 42644328 03/33/7 


SUBROUTINE SRTNOW 


aan ne CURRENT FRAME REPLY PRE=SORT 


iS pra 


COMMON /WINCOW/IWIN, IAZWIN (9) »T WNBEG(9, 3) » THNEND (953) 5 TMNPPT(QD y 
LT MXPRT 
COMMON/SORT/SCRT (2990) »ZSORT (20°09) ySORTNX (2393) ,ISCRTX(2000) ; 
_COMMON/REPLY /REPNOW (1000) NONI JP (1300) yNOWREP,NEXREP, JREPTIM,ISAVE iq 
COMMON/TIME/T IME > TNEXT,FRAM=TyENOTIM yACTINE, IACUPD yFRAMSC 13 
C OMMON/INPUT /COSHBW, REPROB, DPROT y BHOVR2 sNTRACK = }3 
ITRK=6 ia 
00 & M=L,NTRACK 
IF (IAZWIN(M) e£Qe0) GO TO & 
IF (TWNBEG (My 1) oLEsREPTIMeANDoREPTIMeLEsTHNEND(My1)) GO TOS  —_ 
IF (THNBEG (4,2) LESREPTIMSANDoREPTIMeLEeTWNEND(My2)) GO TO S aes 
GO To 4 
ITRK=ITRK +4 
IF (ITRK.GT.5) GO TO 4 
ISAVE=ISAVE*OR*SHIFT (My 364+4*ITRK) 
4 CONTINUE 

IF CITRKe ECO) RETURN 

NOWREP=NOWREP 44 

=(REPTIM@TIME) *FRAMSCH1. 
IF (SORT(K)eLTeO0) GO TO 2 
00 1 M=1,2090 
_KSK+4 

IF (KeGTe2020) K=4 

IF (SORT(K).LT006) GO TO 2 
i C GNTINUE 

PRINT 3 
3 FORMAT (* FRESORT ERROR*) 

sToP 
2 SORT(K)=REFTIM 

ISORT UK) =ISAVE 

RETURN 

£ND 


uu 


nil tlh nt 


LI? 2e FERENCE MAP (R=1) 


>* hE REL CCATION 

3 REAL TIME 3 83WOVR2 REAL INPUT 

b REAL INPUT 2 PROT REAL INPUT 

v 2FAL TIME 2 FRAMET REAL TIME 

c R OAL TIME 5 Iacueo INTEGER TIME 

, INTEGER ARRAY WINDOW 3723 ISAVE INTEGER REPLY 
INTEGER ARRAY SOPT 13560 ISORTX INTEGER ARFAY SCRT 
INTEGER 0 IWIN INTEGER WINDOW 
iNTEGER 61 M INTEGER 

2 ItMTEGER REPLY 175) NOWIJP INTEGER ARRAY REPLY 


F TtTEGER PEPLY 4 NTRACK INTEGER INPUT 


ROUTIN: SRINTX cares ur tne Bese fete ore cen es es 


SUBROUTINE SRTNEX 


a aeee “NEXT FRAME REPLY PRE-SORT _ a 


oq 


COMMCN/WINCOW/IWINy IAZWIN(9) y THNSEG (9, 3)» THNEND (993) y TMNPRT (9) y 
1TMXPRT 
C OMMON/SORT/SCRT (2030) »ISORT (2970) ySORTNX (2900) pISCRTX(20UG) 
C OMMON/RE PLY /FEPNOW (1005) yNOWI J? (1000) yNOW2EP,NEXREFyREPTIM,ISAVE 
COMMON/TIME/ TIME, TNEXT yFRAMET, ENDTIM yg ACTIME y IACUPD yFRAMSC 
: C OMMON/INFLT /COSHBW yREPROBy CPROT y8WOVR2,NTRACK ; P 
ITRK=6 
| 00 4 M=1,NTRACK 
IF (IAZWIN(M) .£Q69) GO TO & 
IF (TWNBEG (1,41) LE@REPTIMSANDeREPTIMSLESTWNENO(My1)) GO TOS 
IF (TWNBEG (M52) 6 LE sREPTIMeAND. O=PTIMSLEeTWNEND(My2)) GO TO 5 : 
TF (THNBEG (443) LE oREPTIMeANDs 2EPTIMeLE«TWN=NO(My3)? GO TO S 
GO TO 4 
5 ITRK=ITRK*+4 
IF (ITRKeGT.5) GO TO & 
CT SAVESISAVE*ORSSHIFT (My 3644* IT PX) 
4 CONTINUE 
IF (ITRK»EG.O) RETURN 
NEXREP=NEXREP +L 
K=(REPTIM*TIME)*FRAMSCHL 
IF (SORTNX(K) LTeGed GO TO 2 
DO 1 M=1,2020 
K=K +4 
IF (KeGT.2039) K=4 
IF (SORTNX(K) LTO.) GO TO 2 
4 CONTINUE 
PRINT 3 
3 FORMAT (* PRESORT ERROR) 
STOP 
2 SORTNX (K) =REPTIM 
I SORTX(K) =ISAVE 
RETURN 
ENO 


SN TYPE RELCCATION 

Powe 2EAL TIME 3 BWOVR2 REAL INPUT 
oa 2 EAL INPUT 2 DPROT REAL INPUT 
ity REAL TIME 2 FRAMET REAL TIME 
vee REAL TIME 5 IACUPD INTEGER TIME 

pi It TEGER ARRAY WINDOW 3723 ISAVE INTEGER REPLY 
T IHTEGER ARRAY SORT 1362 ISORTX INTEGER ARRAY SCRT 

by INTEGER . IWIN INTEGER WINDOW 

| INTEGER 6a M INTEGER 

: INTEGER REPLY 1759 NoWIUP INTEGER ARRAY REPLY 


HUTINE CHRONO TH/ 74 O° T=2 FIN 4.6*4326 03/30/77 


SUBROUTINE CKRONC 
COMMON/SORT/SCRT (29070) sISORT (C2000) ySORTNX (2000) ,ISCRTX(2007) 5 
COMMON/REPLY /REPNOW (1009) ,NOWT UP (1000) yNIWREP,NEXREP,REPTIM, ISAVE 


Tol 


v 
Cress PACK FRE*SORTED ARRAY 
Cc 
K=6 
_.00 2. I=1y,NCWREP_ : : = PAR tee, 
1 K=K+L 


IF (SORT(K)sLT.d.) GO TO 14 
REPNOW (I) =SORT(K) 
NOWIUP (I) =ISORT(K) 
CONTINUE 
Cee eee FINAL BUBBLE SORT a) ee 
M=NOWREP <1. 
IF (MeLEeG) RETURN 
00 6 L=15M 
K=NOWREP@L 
SWAP="1. 
DO 5 I=1,K 
IF (REPNOW(T) -LE*REPNOW(I+1)) GO TO 5 
SWAP=REPNOK(I) 
REPNOW (I) =REPNOW (I+1} 
REPNOW(I#1)=SWAP 
I SWAP=NOWIUP (I) 
NOWIUP (1) =hOWIUP (141) 
NOWIUP (1#4)=1ISWAP 
5 CONTINUE 
IF (SWAP) 40,656 
6 _ CONTINUE 
it RETURN 
ENO 


CLIC REFERENCE MAP (R=1) 


S¥ TYPE RELCCATION 


INTEGER 3723 ISAVE INTEGER REPLY 
T INTEGER ARRAY SORT 13563 ISORTX INTEGER APRAY SCRT 
F INTEGER Sa (K INTEGER 
INTEGER 36 OM INTEGER 
fc INTEGER REPLY 1759 NOWIUP INTEGER ARRAY PEPLY 
ec INTEGER REPLY 0 REPNOW PEAL ARRAY REPLY 
iv ZEAL REPLY 4 §oRT REAL ARKAY SCRT 
hy REAL ARRAY SORT 40 SWAP REAL 
£7°U¢ 
G 32 26 «5 
2t £0 
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UNCLASSIFIED 


AEEEEEE 


ls zg Zs 


Ilzz 


Wes 
lle 


lee 


MICROCOPY RESOLUTION TEST CHART 
NATIONAL BUREAU OF STANDARDS-1963-A 


pQuTImM® LCL2CAR T4s7h OPT=2 FIN 4.644338 63/30, 


SUBROUTINE LL2CAR (RLAT,RLON,H,CAR) 
rn 


uv 


CARnee THIS SUBROUTINE CONVERTS A POSITION SPECIFIED AS LATITUDE, 
GesSs> LEONGITUDE, AND HEIGHT ABOVE A SPHERICAL EARTH TC EARTH=CENTERED 
Estes CARTESIAN COORDINATES (X*AXIS ALONG INTERSECTION OF GREENWICH 
CARES MERIOIJAN WITH EQUATORIAL PLANE, Z*=AXIS ALONG NORTH PCLAR AXIS, 
CAewee YeAXIS CCMPLETES RIGHT@HANDED TRIAD) - 

C 


DIMENSION CAR (3) 
DATA RE434274747/ 
R=RE+H 
RCLAT=R*COS (RLAT) 
CAR(1) =RCLAT*COS (RLON) 
wun MMGUR VANCE ROMS TNCRCONN me 
CAR(3) =R*SINCRLAT) 
RETURN | aad ©: ne NT es Sat ee a 
END 


“AACLIC REFERENCE MAP (R=1) 


NTS i 
ROAR okies a fa BOS Den 2k * 
St TYPE RELCCATION 
R REAL ARRAY FeP. 9 4H REAL FePe 
REAL 27 =RCLAT REAL 
REAL 0 RLAT REAL FePeo 
DIN REAL ; FePe e a te : 
TYPE ARGS 
'S REAL 1 LIBRARY SIN REAL 1 LIBRARY 


AONE RC 


T4/74 oPT=2 


c 
Creve 
Crees 


Cr». te 


a 
y 


INPUT ANGLE ANGe 


SUBROUTINE RCTX (AIN,ANS,AQUT) 


DIMENSION AIN (3) ,AQUT (3) 


__SANG=SIN(ARG) 
CANG=COS (AKG) 
AOUT(4)=AINCL) 
AOUT(2)=AIN(Z)*CANG + 
AQUT(3)=AINC3)*CANG = 
RETURN 
ENO 


T? REFERENCE MAP (R=1)_ 


SN 


bTH 


TYPE RELCCATION 
REAL ARRAY Fes 
REAL ARRAY FePe 
REAL 

TYPE ARGS 
REAL 1 LI@RARY 
238 19 


AIN(3)*SANG 
AIN(2)*SANG 


A-91 


FIN 4644328 


63/30/78 


THIS SUBROUTINE COMPUTES THE TRANSFORMATION ACUT OF THE INPUT 
COORDINATES AIN BY ROTATING POSITIVELY ABOUT TFE X-AXIS 


THROUGH THE 


1 LIBRARY 


| 
| 
| 
| 
| 
| 
] 
| 


ih ins tic eh 


UsnUUILeis KRULL (4/74 OPT=2 FIN 46644328 | 93/3 


' SUBROUTINE ROTZ CAIN, ANG, AOUT) 
G.. ieee een tates: Zs Lo An aee OS Be Aas nO nate Em 
Se eeee THIS SUBROUTINE COMPUTES THE TRANSFORMATION ACUT OF THE INPUT 
CF eR ee COORCINATES AIN BY ROTATING POSITIVELY ABOUT THE Z-AXIS THROUGH 
CHeree INPUT ANGLE ANGe 


© 
¥ 


DIMENSION AIN(3) yAQUT(3) 
_SANG=SINCANG) 
CANG=COS (ANG) 
AOUT(AD=AIN(L)*CANG +# AIN(2)*SANG 
AOUT(2)=AIN(2)*CANG = AIN(1)*SANG 
ACUT(3)=AIN(3) 
RETURN 

END 


PVAOLIC REFERENCE MAP (R=1) 


tvTs 
aTZ 

SN TYPE RELCCATION = Pa 
{N REAL ARRAY SC See se eT TSAR Pee 
qT REAL ARRAY FeP. 22 CANG REAL 
{1G SAL 

TYPE ARGS 

S REAL 1 LIBRARY SIN REAL 1 LIBRARY 
5 3 . : x Gs 


LEtGoty 236 19 


CUTTNE PAGE 74/74 OPT=2 42644338 


03/3577 


SUBROUTINE PAGE 
C Eee BEES bse : 
Ceesee GO TO TOP OF NEXT PAGE 
c 


PRINT 4 
FORMAT (1H1) 
RETURN 


CTINE CUTPUT 74/74 oPT=2 FIN 4e64433E C3/3ur7% 


SUSROUTINE OUTPUT 
COMMON /SOFI/SCRT (2250), 1SORT (297) >SORTNX (2960) gISCPTX(2003) 
COMMON /RE PLY /REPNOW (16Gu) yNOWI JU? (1090) yNOWREP,NEXREF,REPTIN, ISAVE 
COMMON /CONST/ FT2NM y 2EAQTH, THOS 3E 1, HALFPL, THOPIyCTR,RTO Cy ONZOVC 
COMMON/AZ/AZINT 
COMMON SINT ER/NINTS, RHOLIX (79), RHOLTY (79), INTIO (70) »ROTRATI7e) > 
LPRPC7L 9B) gINTRVL (70) gLINT(73),AZSTRT (79) ISTAG (70) 

DIMENSION KTRACK(5) 

DATA MASKS /0009999000030020909178/ 
DATA MASK19/. 5090026600906 19177787 
DATA MASK20/0G900C09C00037777778/ 

A ZOUT=AMOO(AZINT*RTD, 360.) 

PRINT S,NOWREP,AZOUT 

FORMAT (1H#,1I5,* REPLIES, AZ =*F6.1) 
IF (NOWREP*EC.9) RETURN 

DO 2 I=1,NCWREP 

I OUT=SHIFT (NCRIUP (I) »-30) «ANDeMASK10 
JOUT=SHIFT(NCHIUP (I) 5-20) sANO.MASK10 
PCW==.01*FLOAT(NOWIJP(T) «ANC eMASK20) 

_ PRINT LSREFNCHIT) p ENTIRDIOUTHRIOUT 4PO 
FORMAT (1H 9F1i066,14,15,F8.2) 2 
WRITE (10) REPNOW(I) sNOWIJP(T) 
CONTINUE 
RETURN 
ENO 


{TC REFERENCE MAP (R=1) 


REL CCATION 
. AZOUT REAL 
ARRAY oa 7 ¢ REAL. 

FT 2NM REAL 
I INTEGER 
INTEGER : INTER INTRVL INTEGER 
INTEGER ISAVE INTEGER 
INTEGER SORT : IS ORTX INTEGER 
INTEGER ~ INTER JOUT INTEGER 
INTEGER LINT INTEGER 
INTEGER d MASK20 INTEGER 
INTEGER NEXREP INTcGeR 
INTEGER INTER > NOWIUP INTEGER 

REPLY ONEOVC REAL 

COMST D POW REAL 

INTER REARTH REAL 

REPLY 2 RSPTI" REAL 

INTER RHC1 TY SEAL 

INTER RTO REAL 

SORT e SORTNX REAL 

CONST TWO4 Se REAL 


I tT! LAT LONG SCAN AZW) PANG (2) PULSE KE 
1 72 336910 4118 8 15 13-0 342.2 0301450 415-6 1 Qe ( 
2 76 33 55 S7_ 118 24% 23 13¢0 _—«s# OF ot 0801966 536-6 522. 3166 Suet 
3 ¢2 3% 5 45 117 14 12 15.9 203 201134 275. De Ue f 
4 4c wh £2 25 Pte 24 4s 15.0 110.0 0306689 375.6 a6 Ce ( 
5 Se 73 39 4A 117 42 43 13.0 137.8 e201723 390. De 0. i 
6 79 SS GS & “SET £5 36 15-0 299.5 e306 332 2226 Be 0. t 
7 $5 74 F145 117 35 41 13.0 35.5 e301457 300. Je de ( 
& 48s. SF. RG. 4S Ste 20 Pes 600 ___ 22304 so 000691 3705- 92 O. 4 
q 144 St a2 °° £47 55 6 6.0 356-5 e060299 275. Je Oe ( 
10 21 3248 Oo #412 55 e@ 15.0 2498 0002448 275. ae Qe ( 
14 27 3 3 38 413 9 29 5-0 7604 2002336 358. Qe ae 
12 32 32 39 . 411435 6 6.5 310.5 o262327 429. Qe eo ‘ 
13 34 32 39 ¢ #11435 9 660 194.7 0302047 268. Qe Ge t 
46 63 OF tC 0 ER ed ONSET 2... Se He ff 
45 46 35 34° & @2t 4 f 5.0 338.9 0001913 241. de Oe ( 
16 47 a? 22.37 “t27 33 3° 4520 219-6 ~=—§ ee DECT2EF 41d. de Oe 
17 48 37 22 34 129 33 3 15.0 31902 2202258 309-6 Je Ce (i 
18 49 37 27 © 120 34 O 6.0 14267 2001242 25%. 0. Ge ( 
19 S¢ 35-40 3 #2736 ¢ 1523 33502 2962180 395. De Ne { 
20 $1. 9o Gt ¢€ t97 40 0 49.0" 7565 0301489 45G5e De Ge 
24 52 35 46 3 £47 40 15.0 16701 2802220 305. Je Oe r 
22 SE 37 GF. 5 “Set Se. oO 603 169.8 0361570 226.6 de 0. i 
22 cu 33 46 0 “£27 22 © 666 8406 2201145 267. 0. Oe ( 
24 ES SS -47% 3 fi 2a 6 6.9 212.8 09C1959 309. Je Ge c 
25 56 33°18 SS 117 246 90 6.0 178.2 2261110 267. a. a. ( 
26 57 34 S222 417 S4 36 22.0_ 75.6 eUC2305 3805. Be De 
27 53 34 8 & “$27 27 28° “D5s0 556C 2902080 2506.  d« Ue i 
28 €t 76 4E SL 119 43 6 1366 24406 0202116 429. Ze de t 
29 62 34 36 16 447 22 f2 15.5 130.2 eS23145 2676 0-~ 0. t 
36 62 $455 € “StF 23 28 15.0 101.2 e00L6%2 275. 0. Qe ( 
31 €7 37 4S 33 122 7 30 20.6 33.7 e9CC946 30a. 9. 0- { 
32 69 a0 4G 8 ‘fae ft @ 2063 523908 == oe DDGTA2 30050 oe De =k 
32 73 ge 33 S2 227 6 4S 2604 183.7 09041626 250. de 0. ( 
34 72 36 2. 4& 119 54 138 12.0 13901 0396459 3602. 7's 0. i 
35 74 33°47 8c “SSE 3° 6 2U06 396-6 e00C290 300. De Oe ( 
26 75 35 47 § 426 3 6 29.0 260-0 0901775 203. Qe Oe { 
37 77 os &7 fo “fis 24 a 13.9 507 0001603 4d5. De Oe ( 
33 @ 35 5° 32. 118 26 42 oO BHSee = OCG179 276. Oe Oe ft 
39 65 “32 SS 21 4ir 6 45" f3e8 139.5 e3C1057 355. De Ge 4% 
40 86 75 35 16 42150 9 15-66 322.9 0301511 395. Qe O6. { 
ut 87 36 325 52 4121 52 25 15.0 20704 0364191 390.6 de Oe { 
42 Gia 32 52 33 1146 24 S51 5.0 20660 0302270 241. Oe 0. { 
43 eg 37 25 2R 4122 0 56 15.0 28204 0200304 375. Ue 0. ( 
4a ou Zs 4115 118 55 53  G6eC 10004 e201473 296. Ue O6 i 
45 Ct ea at 4i7 42 ¢ 6.0 Sues ef11132 305. a. 0. : 
46 CF a4 3F 53) 418 6 25 15.0 289-6 2101560 270. de Qs ( 
47 c7 35 23 4& 120 21 12 5.0 18569 eMC 365. Ge Qe i 
43 CR 34 216 41417 47 18 10.0 26665 ©201584 296. de Oe ( 
49 129% 3° 7 5 “S95 7 & 15.6 182.6 e301 755 4365 3530 4h. Si 
59 16% % € G7 "439 7 ig 15-0 202.3 2361871 322. 0. Oe ( 
S41 ye GU ¢€. "288 st @ 16.5 174.8 eI09779 36d. Je Oe i 
52 66 f2°A2 3 “S87 2G 4G 1500 28362 031062190 393. 0. Ce { 
53 £07 ee he 6G PET S&S 10.0 32361 eDCC1INF 2656 de Ce i 
5u 458 $2 42 427 164 ¢ 13.0 307-5 ecC21%6 240. ). Oe ‘ 
es 129 st Zé is 149 3a v 15.0 480.5 0901722 360-6 ee Je ( 
56 142 23 42 117 S50 0 6.0 6664 e3GS502 292. Qe Ue ' 
57 416 S640 3° S 6.0 308.9 eOCC745 295, De De \ 
58 $47 Te 4G “ 4if 3 ry 15-0 78-8 Pe Rare 3 1") 2186 Oe Ge : 
A ee 4 Te 3 46" £25 35 S7 6.2 343.2 efC2032 360. 06 Ge 
53 £23 % G42 4125 36 31 15.6 207-7 0001339 275, Ve Ge 
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APPENDIX B 
BCAS MODE SELECTION 


Over 100 modes of BCAS have been analyzed using the 
Simulation tool described in Appendix E. Candidate modes of 


operation were selected based on the criteria of: 


e Accuracy better than 825 feet 
e ATCRBS/DABS/ATARS/ATC compatibility 


s Ability to work everywhere and under all 
conditions 


The initial selection of candidate modes was then configured 
as a preliminary system concept and is illustrated in Figure 
B-1. It should be noted that not all modes meet the required 
position accuracy of 825 feet but are the best available to 
satisfy the additional compatibility and operational criteria 


without adding special features to BCAS. 


It should be further noted that options for addi- 
tional features such as a directional antenna, synchro-DABS, 


azimuth echo and operation with unmodified ATCRBS sites appear 


in this figure. 


The recommended concept, which uses a directional 
antenna and requires synchro-DABS with azimuth echo whenever 
the DABS all-call lockout bit is set, was formulated to elin- 
inate modes with poor accuracy and to provide operational 


capability in high density airspace. Even in this case, a 


few modes exist, particularly in singular regions, where the 


accuracy requirement cannot be met. The penalty in these re- 
gions for reduced accuracy is a higher false alarm rate and 
missed alarm rate but represents a graceful degradation in 


system performance. 


FOR 


FIGURE B-1 


BCAS MODE SELECTION AND CONTROL LOGIC 
(DETAILED FLOW) 


(SEE DRAWING) 


APPENDIX C 


TRAFFIC MODELS AND AIRCRAFT DENSITIES 


A traffic model was required for the technical analysis 


of the various BCAS modes. The L.A. Basin was used for model- 


ling purposes since it represents both the region of highest 


aircraft density, and the largest area of continuous high 

density. Therefore, the Basin represents the pacing terminal 

area for BCAS design. Several traffic forecast models for the 

L.A. Basin in 1982 and beyond were studied (References C-1 and C-2. 
A model (Reference C-1) based on an actual traffic “snapshot" in 


1972 and projected into 1982 was used as the baseline for compariso 


This "full 1982" model of 743 aircraft in a 120 by 120 


nmi square area centered at the L.A. airport was then reduced by 
roughly half to 382 aircraft by randomly deselecting aircraft to 
produce a "half density" model. This was performed again to 

produce a quarter density model of 177 aircraft, so that interim 
measurements and analyses could be accomplished. Examples of the tr 
traffic distribution for each model are illustrated in Figures 


C-1, C-2, and C-3. The quarter density model 
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of 177 aircraft is not significantly different from the 1975 


beacon-only peak traffic count of 150 aircraft. 


In order to simulate a realistic test environment, 
a "target"' aircraft was chosen (#399) which was located in a 


relatively high density (but not the peak density) area of 


the basin (13 nmi East and 15 nmi South of LAX). A ten 
nautical mile radius circle was then drawn on each model 


representation of the basin to indicate where the BCAS air- 


craft simulation was to be tested. The 10 nmi range was used 


to simulate the tracking problems for a BCAS with a 10 nmi 


reception and tracking range. Figure C-4 shows the traffic 


density count taken by the BCAS in 10 degree increments as it 


is positioned around aircraft #399. The counts are given for 


the full, half, and quarter densities counting all the aircraft 


within 10 nmi of the BCAS aircraft. Obviously local peaks occur 
when the BCAS is at roughly 90° and 330° from #399, and a mini- 


mum at 190° (over San Pedro Bay). 


In a separate density analysis, it was assumed that 


any aircraft in the basin could be a BCAS. Each of 743 air- 
craft were examined to determine the number of targets within 
5, 10, 20 nmi with and without altitude filtering of +3000 ft. 


The results of this analysis are presented in histogram form in 


Figures C-5 through C-10. Figure C-1l is a similar histogram but 


shows the number of aircraft around each of the 743 aircraft 
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in the L.A. Basin within 5 nmi and +2000 ft and all others within 


a 60 second volume, determined by relative speed, and closing. 
Figure C-12 illustrates the maximum number of targets within 
10 nmi of another aircraft for the full, half and quarter 


density models. 


Finally, it is possible to relate the LAX model to 
other major hubs in the United States by using Figure C-13. 
In Figure C-13, estimates of aircraft densities for the top 
60 hubs are presented. Note that LAX for 1985 is given as 
0.095 aircraft/nmi-. This figure is based on an assumption 
of uniform density distribution. From Figures C-5 through C-10, 
however, it is observed that peaking factors of 3 or 4 to l 
occur so that local densities on the order of 0.3 to 0.4 air- 
craft/nmi? exist. This has also been verified by recent flight 
tests in LAX where peak densities of approximately 
0.1 aircraft/nm” were experienced of just transponder-equipped 
aircraft. The 1975 figure of 0.048 aircraft /nmi~ in Figure C-13, 
reduced by 30%, approximates the number of transponder aircraft 
¢.039 aircraft/nmi), so that a peaking factor of 3 yields 


« O:1 aircraft/nmi~ which correlates with the test data. 
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XV *SdV ISVOE INNIS 


HUB 1975 
LAX 0.048 
CHI 0,028 
MIA 0,025 
SFO 0, 022 
nyc 0,020 
DFW 0.020 
WAS 0.017 
DET 0.018 
MSP 0.016 
EWR 0,014 
TAH 0.013 
SEA 0.013 
ATL 0.012 
BOS 0.011 
MKC 0, 009 
PHL 0.009 
DEN 0.009 
CLE 0, 009 
PIT 0,009 
STL 0,007 
MSY 0, 006 
LAS 0,004 
BAL 0,003 


HIGH DENSITY 


OKC 0,029 
SAN 0,027 
TUL 0,026 
RIV 0,025 
SIC 0.025 
PDX 0,024 
IND 0.024 
TPA 0,023 
MKE 0,022 
PHX 0,021 
SAC 0,018 
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MEM 0,016 
BDL 0,015 
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SLC 0,008 
ORF 0,008 
ROC 0,007 
ELP 0.006 
RIC 0,006 
RNO 0.005 
BNA 0.005 
SYR 0,005 
CLT 0.005 
RDU 0.004 


1980 1985 


0, 068 0,095 
6,040 0.056 
0,038 0.050 
0,031 0,043 
0,029 0.040 
0,028 0,039 
0,023 0,033 
0, 023 0, 032 
0,022 0,031 
0,019 0,027 
0,019 0,026 
0.018 0,025 
0,017 0,024 
0,016 0,022 
0,013 0,018 
0,01 0,018 
0.013 0,018 
0,012 0,017 
0,012 0,017 
0,010 0,014 
0, 008 0,011 
0,006 0, 008 
0.004 0,006 


1990 


0,133 
0,079 
0,070 
0,061 
0,057 
0.005 
0,046 
0, 046 
0.044 
0.038 
0,037 
0,036 
0,034 
0.031 
0,026 
0.025 
0.025 
0,024 
0.024 
0.020 
0.016 
0.011 
0.009 


HUBS (50nm radius) 


0.041 0,057 
0.039 0,054 
0.036 0.051 
0,036 0,050 
0.035 0,049 
0,034 0,047 
0,033 0,046 
0,032 0,045 
0,031 0,043 
0,029 0,041 
0,028 0,036 
0,026 0,036 
0,025 0,035 
0.022 0,031 
0,020 0,029 
0,020 0,028 
0,019 0,027 
0,018 0,025 
9,017 0,024 
0,017 0,024 
0,015 0,021 
0,014 0,020 
0,014 0,020 
0,014 0.019 
0,013 0.018 
0,013 0.018 
0,012 0,017 
0,011 0,015 
0,011 0.015 
0,009 0,013 
0,008 0,012 
0,008 0,012 
0,008 0,011 
0,007 0,010 
0,007 0,010 
9,007 9,010 
0,006 0,009 


0,081 
0,076 
0.071 
0,071 
0.069 
0.067 
0,065 
0.063 
0,060 
0,057 
0.051 
0.050 
0.049 
0,044 
0,040 
0,040 
0,038 
0,036 
0,034 
0,034 
0,030 
0,028 
0,028 


0,027 
0,026 
0,025 
0,023 
0,021 
0,021 
0,018 
0,016 
0,016 
0,015 
0.014 
0,014 
0,013 
0,012 


MEDIUM DENSITY HUBS (30nm radius) 


Figure C-13: AIRCRAFT DENSITY AT THE 60 TOP HUBS 
(Aircraft per square nautical mile) 
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1995 


0.128 
0,110 
0,098 
0,085 
0,079 
0,078 
0, 064 
0, 064 
0,061 
0,053 
0,051 
0.049 
0,048 
0,044 
0, 0%6 
0,035 
0,035 
0,034 
0,034 
0,028 
0,022 
0,016 
0,012 


0,113 
0,107 
0,100 
0,099 
0,096 
0,093 
0,091 
0. 088 
0,084 
0,080 
0.072 
0,070 
0, 068 
0,061 
0,056 
0,056 
0,053 
0,050 
0.047 
0,047 
0,042 
0,040 
0,039 


0,038 
0.036 
0,036 
0,033 
0,030 
0,030 
0,026 
0,023 
0,023 
0, 021 
0,020 
0,020 
0,019 
0,017 
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APPENDIX D 


RADAR COVERAGE 


The passive BCAS derives its proximity information from 
ground radar interrogations and airborne transponder replies. 
Therefore passive mode BCAS operation is limited to areas of 
suitable ground radar coverage and geometry. This may be 
contrasted with active BCAS operation which derives all of its 


data from air-to-air interrogation and reply. 


COVERAGE 


Low altitude radar coverage is illustrated in Figures 
D-1 through D-4 for the 88 FAA enroute and 162 FAA terminal 


radars in CONUS. The airport surveillance radars are plotted 


for equivalent heights of 1250 and 5000 foot altitudes (50 and 
100 nm ranges). The air route surveillance radars shows 
coverage for 5000 and 20,000 foot altitudes (100 and 200 nm 
ranges). These figures illustrate the vast areas of CONUS 
which have no low altitude radar coverage at all, including 
significant areas in congested regions. Active ranging 
techniques must be employed in areas with no coverage to provide 


adequate collision protection. 


A radar coverage computer program ("CIRCLE") was 
developed to plot the local coverage limit envelopes for radars 
at the top 60 hubs. The program used a file of the locations 
of 698 CONUS radars active as of 1976. The output is a series 
of 120 plots centered at the principal airport of each of the 
60 major hubs at 4000 ft. and 1000 ft. altitudes and are 
included as Attachment 1 of this Appendix. This included the 
23 high density and 37 medium density hubs classified by revenue 
passenger enplanements, and the 301 within those 60 hubs [2]. 

An example is given in Figure D-5 which illustrates the coverage 
limit at 1000 ft. altitude of the four radars in the vicinity 

of the New Orleans' Moisant airport. As in all the other plots, 
the output assumes a smooth spherical Earth line-of-sight 
limitation on radar reception. Terrain and buildings will 
seriously reduce the envelopes but were too complex an issue 

to be dealt with properly within the constraints of this analysis. 
Multipath problems would also limit the usefulness of the 
signals at certain azimuths and elevations. These effects 

were also considered too complex to be treated here, although 
excellent treatments of both ground-to-air and air-to-air 


multipath are available in recent reports (15, 16). 


In Figure D-5 the numbers on the plot are locations of 


the radars, and 45 nm radius circles are drawn about each 


indicating the reception limit at 1000 ft. altitude. There 
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is also a 60 nm radius range ring drawn about Moisant for 
reference, It can be seen that approaching from the Southwest, 
coverage from radar #272 is not obtained until within 7 nm of 

the airvort which allows about two minutes of tracking time. 

It can be seen that the region in Figure D-5 suitable for 
multisite passive BCAS is limited to a roughly 40 by 30 nm region 
inc: jing all the radar positions and located to the Northeast 

of the airport. Any region beyond this will not provide suitable 
geometry and coverage for passive multisite tracking of ground 
interrogations. Also, these systems are susceptible to rapid 
geometric shift when the tracker is within 5 nm of one of the 
radars used for tracking. The geometry in Figure D-5 was 


considered marginal for multisite passive mode tracking. Some 


other mode would be necessary to ensure definite tracking of 


all possible threats within the low altitude maneuvering 
traffic area vicinity. This could take the form of active 

mode enhancement until the garble situation due to traffic 
growth inhibits further user of the active mode. Other means 
include providing a sufficient amount of the threat geometry 
information to the BCAS user via data link to allow the BCAS 

to determine the geometry passively with only single radar site 


information. 


A similar analysis was accomplished for all 60 hubs at 
both 1000 ft. and 4000 ft. altitudes. The assumption was made 
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that the radar deployment would remain constant (i.e., only 
replacing existing sites with newer equipment, no decommissioning 
of old sites or commmissioning of new ones). Figure D-6 


summarizes these results by the total hubs lost at each altitude, 


and the total airports and types of airports within these hubs. 
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APPENDIX E 


BCAS ACCURACY ANALYSIS 


INTRODUCTION 

A universal tool capable of estimating the accuracy of 
any BCAS mode has been developed. This simulation is based on 
"covariance analysis" as opposed to "Monte-Carlo analysis" and 
as such utilizes much less computer time. This simulation was 
used to perform accuracy analysis on over 100 possible BCAS 
modes and to investigate both target and radar singularities 
which occur in certain modes of operation. The results of 
these analyses were used to assist in the selection of the 


recommended BCAS concept. 


DESCRIPTION OF THE SIMULATION 

The mathematical approach used to develop the simula- 
tion tool is based on linearization of the basic equations 
required to solve for target range and bearing relative to 
BCAS. The basic equations which need to be solved are given 


by 


oy" Oye * Pgg ~ Pag 


= Pig cosha, - Pop cos (B-a, 9) 


10 
Pip sinda, = Pom sin (B-a, 9) 


and the geometry for which these equations apply is illustrated 


in Figure E-l. 


Note that these equations are valid only for a two- 
dimensional analysis, but the results of this two-dimensional 
analysis have been compared with a complete three-dimensional 
Monte-Carlo simulation of a specific mode of BCAS. This com- 
parison of the two analytic methods showed equivalence of the 
approaches. Also note that a covariance analysis of the three- 
dimensional equations is also possible, but was not done as 
part of this work due to time limitations. This will be pur- 


sued at a later time. 


Linearization of equations (E-1l), (E-2), and (E-3) 


results in the following equation set: 
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0 sindo., Por* Pip 9 -sin (8-4, ) ~Pop* 88 0 
cos (8-a., 4) cosha., cos (B-a, 4) 
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Equation set (E-4) clearly represents three equations in seven 
unknowns. Thus, to obtain a deterministic solution, four 
variables must be specified or measured. An over-determined 


solution occurs when more than four variables are measured. 


The total number of possible modes of BCAS can then 


be calculated from 
N=C, + CL + c! + Cc E-5) 
6 ; ( 
so that N = 64. 


The simulation program "SINGLE" has the capability of 
analyzing any of the 64 possible solutions for any radar/BCAS/ 
target geometry. The technique used is to form the appropri- 


ate matrix equation from (E-4) in the form 
A {ox} = B {dy} (E-6) 


where {6x} are the unknown variables and {sy} are the 
measured variables. The solution of (E-6) in general is given 
by 

éx = [Lata] ~! at™ B oy (E-7) 


The covariance matrix of the unknowns is thus related to the 


covariance matrix of the measurements by 


P=KMK (E-8) 
where 
k = (ata)? a™ B 
M = E (éy éy’) (E-9) 
= E (6x 6x?) 


Given P, the main diagonal elements are the variances of the 


unknowns. 


As an example, suppose a BCAS mode which measures 


10’ ba, and Th Under this assumption, 


5P 19 é 
Sao a 
6y = and x = 50 on 
Sda, 
68 
ét, 


M is then formed by placing the variance of the measurement 


on the main diagonal so that 


Z 
Oo 
Pio : ° ° . 
Oo 
0 X16 0 0 
M = o2 
0 0 Aa, 0 
32 | 
0 0 0 T, 
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K is solved for using specific geometry and then P is calcu- 


lated using equation (E-8). The position uncertainty is then 


obtained from 


o 2 2 2.4 - 
Op = {o a oat (E-10) 


This technique has been extended to two radars and is 
documented in program "DUAL." A third program called 
"DUAL-ATCRBS" has the capability of analyzing two-radar/two- 
target BCAS modes. Listings of these programs are presented 


as Attachment 1 of this appendix. 


E.3 RESULTS 

The basic geometry used for the accuracy analysis is 
illustrated in Figure E-2. The results of these analyses are 
presented as Attachment 2 at the end of this appendix. A plot 
of position uncertainty vs. target bearing and the associated 
numerical data is presented. Range, bearing, and total 
position uncertainty as a function of target bearing are 


contained in the data. 


E.4 CONCLUSIONS 
A clear advantage of this approach is that a specific 
computational algorithm for a given BCAS mode does not have to 


be developed for purposes of accuracy analysis. A Monte-Carlo 
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simulation has been developed which has the ability of 


evaluating a specific BCAS computational algorithm and has 
been used to check the validity of the covariance analysis. 

As expected, both results agree to a high degree of accuracy. 
The listing of the Monte-Carlo simulation is presented as 
Attachment 3 at the end of this appendix. Although this 

work is two-dimensional, it is relatively straightforward to 
extend to three dimensions. It is planned that this work will 
be done in the near future; it will be used to investigate the 
possibility of estimating target altitude when the target is 
equipped only with a Mode A transponder (without an encoding 


altimeter). 
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RECOMMENDED CONCEPT 


ATTACHMENT II 


Mode 


Single # 


Single # 


Single # 


Single # 


Single # 


Dual #13 


Dual #20 


Dual #26 
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9 


15 


ES 


SUMMARY CHART 


Range 


BCAS/Target 


10/5 
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50/10 
100/20 
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50/10 
100/20 
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ATCRBS 


Semi-active 
ATCRBS and DABS 


Semi-active ATCRBS 


Semi-active ATCRBS 
with Directional 
Antenna with 1° Error 


Semi-active 


ATCRBS and DABS 


Semi-active ATCRBS 
Passive 


Passive 


Semi-active 


SUMMARY CHART (CONTINUED) 


Range 
Mode BCAS/Target Comments 
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